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ABSTRACT 

Quantitative description of the statistics of intensity fluctuations within spec- 
tral line data cubes introduced in our earlier work is extended to the absorbing 
media. A possibility of extracting 3D velocity and density statistics from both 
integrated line intensity as well as from the individual channel maps is analyzed. 
We find that absorption enables the velocity effects to be seen even if the spectral 
line is integrated over frequencies. This regime that is frequently employed in 
observations is characterized by a non-trivial relation between the spectral index 
of velocities and the spectral index of intensity fluctuations. For instance when 
density is dominated by fluctuations at large scales, i.e. when correlations scale as 
r -7 , 7 < 0, the intensity fluctuations exhibit a universal spectrum of fluctuations 
~ K~ 3 over a range of scales. When small scale fluctuations of density contain 
most of the energy, i.e. when correlations scale as r -7 , 7 > 0, the resulting 
spectrum of the integrated lines depends on the scaling of the underlying density 
and scales as K~ 3+1 . We show that if we take the spectral line slices that are 
sufficiently thin we recover our earlier results for thin slice data without absorp- 
tion. As the result we extend the Velocity Channel Analysis (VCA) technique to 
optically thick lines enabling studies of turbulence in molecular clouds. In addi- 
tion, the developed mathematical machinery enables a quantitative approach to 
solving other problems that involved statistical description of turbulence within 
emitting and absorbing gas. 
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1. Introduction 

There is little doubt that interstellar medium is turbulent (see reviews by Armstrong, 
Spangler & Rickett 1995; Lazarian 1999a, Lazarian, Pogosyan & Esquivel 2002). Turbulence 
proved to be ubiquitous in molecular clouds (Dickman 1985), diffuse ionized (Cordes 1999) 
and neutral (Lazarian 1999b) media. Magneto-hydrodynamical (MHD) turbulence controls 
many essential astrophysical processes including star formation, transport and acceleration of 
cosmic rays, of heat and mass (see Schlickeiser 1999, Vazquez-Semadeni et al 2000, Narayan 
& Medvedev 2001, Wolfire et al. 2002, Cho et al. 2003, see recent reviews by Cho, Lazarian 
& Vishniac 2002, henceforth CLV02a, Elmegreen & Scalo 2004, Lazarian & Cho 2004). 
Unfortunately, we are still groping for the basic properties of interstellar turbulence. For 
instance, it is unclear whether turbulence in molecular clouds is a part of a global turbulent 
cascade or has its local origin. 

Recent years were marked by a substantial progress in theoretical and numerical de- 
scription of both incompressible and compressible MHD turbulent cascade (Goldreich & 
Shridhar 1995, Lithwick & Goldreich 2001, Cho & Lazarian 2002a, CLV02a, Cho & Lazar- 
ian 2003ab). It is encouraging that the measured statistics of electron density fluctuations 
(see Armstrong et al. 1995), synchrotron and polarization degree fluctuations (see Cho & 
Lazarian 2002b) are in rough correspondence with the theoretical expectations (see more 
discussion in Cho & Lazarian 2003b). Further testing requires better input data, and this 
makes the determination of actual power spectra of interstellar turbulence extremely impor- 
tant and meaningful. Indeed, it gives a chance to distinguish between different pictures of 
turbulence, e.g. shocks versus eddies picture, and put theoretical constructions developed for 
various interstellar processes, e.g. cosmic ray propagation (see Yan & Lazarian 2002), grain 
dynamics (Lazarian & Yan 2002, Yan & Lazarian 2003), stochastic reconnection (Lazarian & 
Vishniac 1999, Lazarian, Vishniac & Cho 2003) on a solid ground. Even within the picture 
of MHD cascade that is advocated by the simulations in Cho & Lazarian (2002) and which 
does not contain shocks, the distribution of energy between Alfven, slow and fast modes does 
depend on the turbulence driving. Therefore the spectra measured from observations can 
be used to constrain the driving. In addition, the particular features of the power spectrum, 
which are associated with energy injection and dissipation are important. Identification of 
such features and the spatial variations of their properties will shade light to the dynamics 
of interstellar gas and processes of star formation (see Lazarian & Cho 2003). 

We should stress that while the statistics of density fluctuations can be only considered 
as an indirect way of testing turbulence, the spectral surveys contain the information about 
the direct measure, namely, the velocity statistics. This statistics is extremely valuable, 
provided that we can extract it from the observational data. As MHD turbulence is indeed 
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an interdisciplinary subject, high resolution studies of interstellar turbulence can test con- 
temporary ideas of MHD cascade. The implications of the improved understanding would 
affect description of wide range of astrophysical phenomena from solar flares to gamma ray 
bursts (see Lazarian et al. 2003). 

Additional motivation for studying interstellar turbulence stems from the advances in 
direct numerical simulations of interstellar medium (see review by Vazquez-Semadeni et al 
2000). Unlike MHD simulations focused on obtaining fundamental properties of turbulence 
(see CLV02a) those of interstellar medium attempt to include various aspects of interstellar 
physics. This usually limits the inertial range covered. However, testing of the numerical 
results against observations is becoming more interesting as numerical resolution increases. 
With other techniques of comparing numerics with observations having their problems (see 
review by Ostriker 2002, also Brunt et al. 2003) power spectra and second order correlation 
functions look as a promising tool for the future (see Lazarian 1995, Lazarian & Pogosyan 
2000, henceforth LP00, Vestuto, Ostriker & Stone 2003). The potential of combining different 
techniques in order to get the properties of turbulence is also large. 

Although in this paper we primarily refer to interstellar processes, the technique we 
discuss has a broader application. For instance, recent advances in X-ray astronomy have 
allowed to get information on turbulence in intracluster gas (see Inogamov & Sunyaev 2003, 
Sunyaev, Norman & Bryan 2003). The issues that researchers face there are similar to those 
dealt in interstellar turbulence studies. 

To get insight into turbulent processes statistical approach is useful (see Monin & Ya- 
glom 1976, Dickey 1995). So far, in astrophysical context the most tangible progress was 
achieved via scintillations and scattering technique (see Spangler 1999). Those measure- 
ments are limited to probing fluctuations of electron density on rather small scales, namely, 
10 8 -10 15 cm. This research profited a lot from an adequate theoretical understanding of 
processes of scintillations and scattering (Goodman & Narayan 1985) and this made it very 
different from other branches of interstellar turbulence research. In comparison, formation of 
emission line profiles in turbulent media has not been properly described until very recently 
and it is natural that numerous attempts to study turbulence in diffuse interstellar medium, 
HII regions and molecular clouds using emission lines (see Munch 1958, O'Dell 1986, O'Dell 
& Castaneda 1987, Miesh & Bally 1994, reviews by Scalo 1987, Lazarian 1992) were only 
partially successful. This is very unfortunate as the line profiles contain unique information 
about turbulent velocity field. We reiterate that the studies of stochastic density provide only 
indirect insight into turbulence and cannot distinguish between active and fossil turbulence 
pictures. 

Studies of the velocity field have been attempted at different times with velocity cen- 
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troids (e.g. Munch 1958, Miesh, Bally & Scalo 1999, Miville-Deschenes, Levrier & Falgarone 
2003). However, it has long been realized that the centroids are affected, in general, by both 
velocity and density fluctuations (see Stenholm 1989). An important criterion when centroids 
indeed reflect the velocity statistics was obtained in Lazarian & Esquivel (2003, henceforth 
LE03). LE03 showed that when the criterion is not satisfied the centroids are dominated by 
density. LE03 pointed out that the velocity centroids can be modified to extract the velocity 
contribution (LE03), but the resulting measures may be pretty noisy. As the result, velocity 
centroids taken alone may not be sufficiently robust and reliable technique. A recent study in 
Esquivel & Lazarian (2004) revealed that for Mach numbers larger than 2 or 3 the traditional 
velocity centroids fail. As Mach numbers as high as 10 are expected for molecular clouds 
the applicability of the traditional centroids to molecular clouds is highly questionable. 

Principal Component Analysis (PCA) of the emission data (Heyer & Schloerb 1997, 
Brunt & Heyer 2002ab) has been suggested as a new way to study interstellar turbulence. 
However, recent testing showed that it provides statistics different from power spectra (Heyer, 
private communication, Brunt et al 2004). What sort of statistics we can obtain using the 
PCA requires a further study. 

Wavelet analysis (see Gill & Henriksen 1990), Spectral Correlation Functions 1 (see 
Rosolowsky et al. 1999, Padoan, Rosolowsky & Goodman 2001), genus analysis (see Lazar- 
ian, Pogosyan & Esquivel 2002) are other important statistical tools. They can provide 
statistics and insight complementary to power spectra (see review by Lazarian 1999a). Syn- 
ergy of different techniques should provide the necessary insight into turbulence and enable 
the comparison of observational statistics with theoretical expectations (see Cho & Lazarian 
2003a for a review). 

Velocity statistics contributes to intensity variations observed in channel maps. Power- 
law spectra suggestive of underlying turbulence were obtained at different times by a number 
of researchers (see Kalberla & Mebold 1983, Green 1993). The problem that plagued those 
studies was inability to separate contributions due to velocity and density. Indeed, both 
velocity and density fluctuations affect the small-scale emissivity fluctuations observed at a 
given velocity. Therefore the separation of contributions to velocity and density requires a 
quantitative description of spectral line data cube statistics. 

This problem was addressed in Lazarian & Pogosyan (2000, henthforth LP00), where 



1 Spectral Correlation Functions if they are measured for the same velocity present just another way of 
dealing with fluctuations within channel maps. Therefore all the VCA results are applicable to them. If one 
tries to generalize Spectral Correlation Functions to the velocity direction (see Lazarian 1999) only limited 
information about turbulence can be available using the tool (see Appendix C). 
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the fluctuations of intensity in channel maps were related to the statistics of velocity and 
density. LPOO introduced a new statistical technique that was termed in Lazarian, Pogosyan 
& Esquivel (2002) Velocity Channel Analysis (henceforth VCA) . Within VCA the separation 
of velocity and density contributions is obtained by changing the thickness of the analyzed 
slice of the Position-Position- Velocity (PPV) data cube. The theoretical description of the 
intensity statistics for both thick and thin velocity slices enabled LPOO to disentangle velocity 
and density contributions to intensity fluctuations in PPV cubes. 

The VCA was successfully tested numerically in Lazarian et al. (2001) and Esquivel 
et al. (2003). Those tests employed the velocity and density obtained via simulations of 
compressible MHD turbulence to obtain synthetic maps, which were analyzed using VCA to 
recover underlying velocity statistics. Since then the VCA has been successfully applied to 
get velocity spectra of turbulence (see section 5). 

The VCA in its existing form, however, does not account for absorption of radiation. 
LPOO argued that for the 21 cm emission data the effect of absorption is marginal in the 
direction of the galactic anti-center for which the technique was applied. However, the data 
in Dickey et al (2001) show substantial absorption in HI towards inner part of the Galaxy. 
Moreover, the absorption is definitely essential for other, e.g. CO transitions. 

The purpose of this paper is to extend the earlier obtained theoretical description of 
the Doppler-shifted spectral line data cubes by including the treatment of absorption effects 
and thus to obtain a quantitative tool to study turbulence in various interstellar conditions, 
including molecular clouds. We introduce the 3D anisotropic statistics of fluctuations within 
PPV data cubes in section 2, describe effects absorption in section 3, provide the statistics for 
thin and thick slices in section 4, explain some of our findings in physical terms in section 5, 
discuss our results and their implications for observational data in section 6. The summary 
is provided in section 7. The lengthy derivations are collected in the appendixes which are 
an important part of the paper. 



2. Statistics of Velocity and Density 

2.1. Basics of Turbulence Statistics 

Statistics in Real Space 
In the presence of magnetic field MHD turbulence gets axisymmetric in the system of refer- 
ence related to the local direction of magnetic field. In this system of reference magnetic fields 
can be easily mixed by turbulence in the direction perpendicular to magnetic fields. Those 
mixing motions generate wave-like perturbations propagating along magnetic field lines. As 
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the kinetic energy decreases with the decrease of the scale, the weak mixing motions bend 
magnetic field lines less and less and the eddies get more and more elongated. 

The Goldreich & Shridhar (1995) model of incompressible turbulence prescribes Kol- 
mogorov scaling of mixing motions perpendicular to magnetic field lines (see a simplified 
discussion in CLV02a) and the elongation of the velocity fluctuations that increases as Z 1 / 3 . 
Further research in Lithwick & Goldreich (2001) and Cho & Lazarian (2002) has shown that 
the basic features of the Goldreich-Shridhar turbulence carry on for Alfvenic perturbations 
to the compressible regime 2 . 

Observations, however, are usually unable to identify the local orientation of magnetic 
field and deal with the magnetic field projection integrated over the line of sight. As the result 
the locally defined perpendicular and parallel directions are mixed together in the process of 
observations (see CLV02a). The fluctuations with more power dominate the signal for both 
the direction of the averaged B and perpendicular to it. There is some residual anisotropy, 
but this anisotropy is scale-independent and is determined by the rate of the meandering 
of the large scale field. In other words, from the observational point of view the spectra of 
intensity fluctuations obtained from the Goldreich & Shridhar turbulence and the isotropic 
turbulence are similar. Observations in Green (1993) and numerical studies in Esquivel et 
al. (2003) support this picture. 

The facts above allow for a substantial simplification of the turbulence description, 
namely, they permit us to use standard isotropic statistics (Monin & Yaglom 1976). Ob- 
taining this statistics would correspond to obtaining the scaling of density and velocity 
fluctuations perpendicular to magnetic field lines. 

Statistically isotropic in xyz-spa.ce density field p(x) has the correlation function 

£(r)=£(r) = (p(x)p(x + r)> . (1) 
We shall also use the correlation function of density fluctuations 5p = p — (p) 

f (r) = (5p(x)Mx + r)> = £(r) - (p) 2 . (2) 

The structure function 

rf(r) = ((p(x + r)-p(x)) 2 ) , (3) 



2 Magnetosonic fast modes that arise in compressible fluid, are, however, isotropic (Cho & Lazarian 2002, 
2003ab). This entails many important astrophysical consequences (see review by Lazarian, Cho & Yan 2002). 
When both fast and Alfven modes are present in the fluid, the total anisotropy of magnetic and velocity 
fluctuations decreases. 
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is another way of describing turbulence. Here d(r) for fluctuations coincides with d(r). 

Statistical descriptors of fluctuations can often be assumed to have power-law depen- 
dence on scale oc r~ 7 . If the correlation is dominated by large scales, the structure functions 
are used and 7 < 0, while for small scale dominated statistics the correlation functions 
are usually used and 7 > 0. We show in Appendix A that the two cases in the asymp- 
totic regime of small r can be treated very similarly. Therefore we will frequently refer to 
correlation functions having in mind both 7 > and 7 < 0. 

An isotropic velocity field u(x) is fully described by the structure tensor (AiijAuj), 
which can be expressed via longitudinal Dll and transverse D^n components (Monin & 
Yaglom 1972) 

(Au t Auj) = {D LL (r)-D NN {r)) r -^ + D NN {r)5 tk , (4) 

where 5^ equals 1 for i — k and zero otherwise. We define z-projection of the velocity 
structure function as 

D z (r) = (AuiAuj)ziZj = D NN (r) + [D LL (r) - D NN (r)] cos 2 9 , cos6> = f • z (5) 
which for the power-law velocity gives 

Tfl 

D z (r) = Cr m [l + -(l-cos 2 9)} , (6) 

if we assume that the velocity is solenoidal 3 . For this paper it is only important that 
D z ~ Cr m . A discussion of the shallow and steep spectra of density as well as of the 
velocity statistics within finite size clouds is given in Appendix A. 

Spectra and Correlation Functions 
Spectra and correlation/structure functions are two complementary ways of describing tur- 
bulence. Given an N-dimensional correlation function £jv( r ) one can obtain the spectrum 

P(k) = J ^re lkr 6v(r) , (7) 

where the integration is performed in the iV-dimensional space. 

In LP00 the Fourier space statistics, namely, spectra were widely used with the cor- 
relation functions playing an auxiliary role. In the present paper we deal with absorption 



3 Separating solenoidal and potential components of the velocity is an important problem that we do not 
address here. It was suggested in LE03 that this separation could be done by combining VCA and velocity 
centroids. 
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defined in real space. Therefore we use real space statistics, namely, correlation and struc- 
ture functions functions. It is obvious from eq. (7) that for power-law correlation function 
(CF), i.e. £at ~ the spectral index is also a power-law, i.e. P ~ k n , where n = —N + 7. 
In other words, 

(spectral index) = (—dimensions of space — CF index) (8) 

In Kolmogorov turbulence passive scalar density correlations scale as r 2//3 , which corresponds 
to 7 = —2/3 in our notations. Thus the Kolmogorov spectrum index is —11/3. In turbulence 
literature E(k) = 47rk 2 P(k) is usually used. In this notation the Kolmogorov spectrum is 
E(k) ~ /c~ 5 / 3 , and thus, the Kolmogorov law is often referred to as —5/3 law. 



2.2. Statistics in PPV 

One does not observe the gas distribution in the real space galactic coordinates xyz 
where the 3D vector x is defined. Rather, intensity of the emission in a given spectral line 
is defined in Position-Position- Velocity (PPV) cubes towards some direction on the sky and 
at a given line-of-sight velocity v . 4 In the plane parallel approximation the direction on 
the sky is identified with xy plane where the 2D spatial vector X is defined, so that the 
coordinates of PPV cubes available through observations are (X, v) . The relation between 
the real space and PPV descriptions is defined by a map (X, z) — > (X, v). 

The central object for our study is a turbulent cloud in PPV coordinates. The statistical 
properties of the PPV density p s (X, v) depend on the density of gas in real galactic coordi- 
nates, but also on velocity distribution of gas particles. Henthforth we use the subscript s 
to distinguish the quantities in (X, v) coordinates from those (X, z) coordinates. We shall 
always assume 2D statistical homogeneity and isotropy of p s (X, v) in X-direction over the 
image of a cloud. However, homogeneity along the velocity direction can only be assumed 
after additional considerations, if at all. Naturally, there is no symmetry between v and X. 

We shall use the following notation for statistical descriptors of the density in PPV 
space: the mean density is p s (vi) = (p s (Xi, Vi)), and the correlation functions of the density 
and, closely related, the density fluctuations 5p s (X 1 , v\) = p s (X 1 ,vi) — p s ( v i) are 

t s {R,v 1 ,v 2 ) = (p s (X 1 ,v 1 )p s (X2,v 2 )) (9) 
is(R,v 1 ,v 2 ) = (Sp s (X 1 ,v 1 )5p s (X 2 ,v 2 )) = $,(R,v 1 ,v 2 ) - p s (vi)p s (v 2 ) . (10) 



All velocities are the line-of-sight velocities, we omit any special notation to denote z-component. 
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Here we have maintained notation which highlights the symmetries of the statistics. Homo- 
geneity and isotropy in Position-Position coordinates X lead the mean density to depend 
only on velocity while the correlation function depends only on the magnitude of separation 
between two sky directions R = |R| = |X X — X 2 |. Dependence on the velocity is retained 
for the time being in the general form. 

For PPV statistics it is more convenient to use the structure functions 

d s (R,v 1 ,v 2 ) = ([p a (Xi,ui) - p s (X 2 ,v 2 )] 2 } (11) 
d s (R,v 1 ,v 2 ) = ([5p s (Xi,t>i) - 5p s (X 2 , v 2 )f) = d s (R, v u v 2 ) - \p s {vi) - p s (v 2 )} 2 (12) 



Let us derive the two-point correlation function in PPV space without using the power 
spectrum formalism employed in LPOO. The line-of-sight component of velocity v at the 
position x is a sum of the regular gas flow (e.g., due to galactic rotation) v ga i(x), the turbulent 
velocity u(x) and the residual component due to thermal motions. This residual thermal 
velocity v — f 9a z( x ) ~~ M ( x ) has a Maxwellian distribution 

(v -v ga i(x) -«(x)) 2 " 



(x)df 



(2tt/5) 1 /2 



cxp 



2(3 



dv 



(13) 



where (5 = kbT /m, m being the mass of atoms. The temperature T can, in general, vary 
from point to point. 



The density of the gas in the PPV space p s (X, v) can then be written as 



p s (X, v )dXdv 



dz p(x)0„(x) 



(14) 



where p(x) is the (random) density of gas in real space. This expression just counts the 
number of atoms along the line-of-sight that have z-component of velocity in the interval 
[v, v + dv]. The limits of integration are defined by the spatial extend S of the emitting gas 
distribution. 

To compute the correlation function £ s we disregard the effect of correlations between 
density and velocity fluctuations 5 : 

(p s (X 1 ,v 1 )p s (X 2 ,v 2 )) = dz! dz 2 (p(xi)p(x 2 )) (0 v i(xi)0 w2 (x 2 )) , (15) 

Jo Jo 

(p s (X,v)) = [ ^<p(x))<^(x)> (16) 

^0 



5 A similar assumption was used in LPOO. It was tested by Lazarian et al (2001) and Esquivel ct al. 
(2003), that this assumption does not cause any problem velocity and density are correlated. 
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The brackets (• • •) in eq. (16) denote statistical averaging over realizations of the random 
density p(x) and the turbulent velocity tt(x) of gas. In particular, for the density we can 
write (p(x)) = p and (p(xi)p(x 2 )) = £(r), without any assumptions on density statistics. For 
the Gaussian velocity field u, described by structure functions given by eq. (4), the statistical 
average of the Maxwell functionals (13) is given in Appendix B (see (B8) (B9), (B4)). 

The general expression for the correlation function given in the Appendix B (see eq.(B3) 
can be simplified if one neglects the edge effects associated with the wings 6 of spectral lines 
(see derivation of eq. (B7)). In this case we have 



2n(CS)^ j_ s j lzl/2 



dz 



S-\z\/2 



dz + £(r) [ J D,(r) + 2/?]- 1 / 2 exp 



(i 



J gal, 



2{D z {v) + 2(3) 



;i7) 

where eq. (A2) was used to write the prefactor explicitly. For a homogeneous statistics all 
correlation functions depend only on the separation between points. Here and further on 
we use variables without indexes to denote separation, as in R = X x — X 2 , z = z 1 — z 2 , 
v — v± — V2 and v ga i = v ga i t i — v ga i^- Variables with an index '+' denote the average of 
two arguments z + = {z\ + z 2 )/2, etc. Dependence on these last quantities appears only if 
statistical homogeneity is broken, e.g. by a finite spatial extend of the gas distribution. 

In the absence of the turbulent and thermal motions, D z ,j3 — > 0, the kernel in eq. (17) 
reduces to the delta-function 5{y — v ga i), and for monotonic v ga i(z) one can fully recover the 
line of sight position of an emitting atom from its velocity. In a general case, while thermal 
broadening just smothers fluctuations out, the effect of the turbulent velocity fluctuations is 
scale dependent and leads to a change of the correlation function slope. 

Equation (17) allows for an arbitrary form of the regular flow v ga i- A variable regular 
flow produces a complex map from galactic to velocity space by itself. Two tractable cases 
are when the regular flow is of a simple linear form, or the one when it is possible to attribute 
all motions to turbulence, thus setting v ga \ = 0. In LP00 we have considered the case of 
linearized shear flow v ga i = f~ 1 z due to Galactic rotation: 



UR,v) ~ / dz f (r) [D z (r) + 2/3]" 1 / 2 exp 



(v-f- l z) 2 
'2(D z (r) + 2/3). 



(18) 



where we have omitted an unimportant dimensional prefactor (see eq. (B9) for the explicit 
form of the prefactor). Galactic shear introduces a natural scale A = [pC]* 11 ™ , at which the 
velocity dispersion ~ C\ m becomes equal to the squared difference of the regular velocities 



6 Edge effects make the image of turbulence inhomogeneous. It is intuitively clear that those effects should 
not affect the statistics on the scale much less than the spectral line width. 
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determined by Galactic rotation (i.e. f~ 2 \ 2 ). Asymptotics obtained in LPOO correspond 
to the scales I which are much smaller than A. Although at these scales galactic rotational 
velocity is much larger than the turbulent one, it's gradient is much smaller than the tur- 
bulent one. It is because of that the LPOO results in the presence of Galactic rotation and 
without it coincide. The regime when the small scales turbulent shear exceeds the regular 
one is natural. Indeed, flows with high Reynolds number produce turbulence with shear 
larger that that of the original flow. This is true, for instance, for Couette flows in the pres- 
ence of the moderately strong magnetic field (Velikhov 1959, Chandrasechar 1962, Balbus & 
Hawley 1991). The corresponding criterion on the magnetic field is usually satisfied in the 
interstellar medium. 

If gas is confined in an isolated cloud of size S and the galactic shear over this scale is 
neglected we get 



v 



UR,v) ~ J* dz(l - M) rn J (r Li/ a ex P ' 



2(£> z (r) + 2/3)J 



(19) 



. 5 v 2S J [£>,(r) + 2/3]V2 

where again somewhat ugly prefactor entering eq. (B9) is omitted. Turbulent effects remain 
important up to the scale of the cloud S at which turbulent structure function saturates at 
the value ~ (CS m ) 1 / 2 . The cloud size S now plays the role of the scale A. This observation 
allows to translate the LPOO results mostly written for Galactic HI with rotation curve 
mapping to a case of an individual cloud. Rigorous calculations provided in LPOO prove that 
for scales much smaller than S one can use the results obtained for infinite medium in the 
presence of shear and substitute A for S. 

When the amplitude of fluctuations grows with separation, one should use structure 
functions in PPV. The transfer from one type of statistics to another is similar to that 
discussed in Appendix A. 

Important feature of the PPV space is that p s (X, v) exhibits fluctuations even if the 
flow is incompressible and no density fluctuations are present. Indeed, when one substitutes 
the expanded expression (A4) £(r) = p 2 + p 2 (r /r) 7 into eqs (17-19), both terms will give rise 
to non-trivial contributions to £ S (R, v). We shall therefore split the result correspondingly 

l(R,v) = i v (R,v)+i p (R,v), 

d s (R,v) = d v (R,v)+d p (R,v). (20) 

with the f-term describing pure velocity effects, while the p-term arising from the actual 
real space density inhomogeneities that are modified by velocity mapping. To simplify the 
notation we have dropped index s from the right-hand-side quantities, since this split is only 
meaningful in PPV space. In Appendix C we discuss asymptotic small-R scalings of different 
contributions to PPV structure and correlation functions. 
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3. Turbulent Statistics and Radiative Transfer 

We start with the standard equation of radiative transfer (Spitzer 1978) 

dl v = -g v I v ds + j u ds . (21) 

In the case of self-adsorbing emission in spectral lines that is proportional to first power of 
density: 

g v = a(z)p(z)(f) v (z) , (22) 
3v = £p( z )M z ) , ( 23 ) 

where (f> v (z) is given by (13). 

A solution of this equation if no external illumination is present is 

I v = e dzp(z)</> v (z)exp(- a(z')p(z')(/> v (z')dz'^ . (24) 

To integrate (24) we shall assume that a is constant. This constancy is the essence of the 
Sobolev approximation that was found to be useful in many astrophysical applications. In 
this approximation we can use an integration variable 

Y v (z) = [ p{z')^ v {z')dz' , (25) 
Jo 

which value at x — S coincides with the density in PPV coordinates, Y V (S) = p s (X, v), and 
Y v (0) = 0, to integrate eq.(24) 

IJX) = e r dY v e- aY " = - \l - e - ap ^' v) ] . (26) 
Jo « 

In the case of vanishing absorption, the intensity is given by the linear term in the expansion 
of the exponent in eq. (26) 

I v (X.) = ep s (X,v) . (27) 

and reflects the PPV density of the emitters. If, however, the absorption is strong, the in- 
tensity of the emission is saturated at the value e/a wherever p s (X, v) ^> 1/a. Identification 
of the low contrast residual fluctuations may be difficult in practice. 

The approximate expression given by eq. (26) can be used to find the statistics of the 
expected emissivity fluctuations. First of all, the mean profile of the line is given by 



(I V (X)) = - [1 - (e-""^))] . (2* 
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Our main goal is to calculate the structure function of the observed emissivity 

P(R) = ^ (/ IvQtx)W v dv - j I v (X 2 )W v dv^j ^, R = Xi - X 2 . (29) 

Here the window function W v describes how the integration over velocities is performed. 
When W v = 1 the integration is being performed over the whole line as it is a frequent case 
for CO turbulence studies (see Falgarone et al. 1998, Stutzki et al. 1998) while measurements 
in velocity slices of PPV data cube (channel maps) correspond to W v strongly peaked at a 
particular velocity. Within the VCA technique by varying the width of velocity channels one 
can obtain statistics of turbulent velocities and density inhomogeneities of emitting medium. 
Note, that the minimal width of the velocity channel is determined by the resolution of an 
instrument. 

Using elementary identities to transfer from averaging over 2D intensities to averaging 
the underlying 3D fluctuations (see Lazarian 1995), in the limit when the absorption can be 
neglected we rewrite the expression (29) in the form 

P(R) = e 2 J dv 1 W(v 1 ) J dv 2 W(v 2 ) [d s (R, Vl , v 2 ) - 4(0, v u v 2 )\ (30) 

where homogeneity in physical directions is assumed. For infinite emitting medium homo- 
geneous turbulence produces a homogeneous image in the velocity space. For finite emit- 
ting cloud the PPV image is approximately homogeneous over velocity separations much 
less than the Doppler linewidth. Notably, the combination in brackets does not depend 
on the mean Doppler broadened profile of the spectral line, d s (R, Vi,v 2 ) — d s (0,vi,v 2 ) = 
d s (H,Vi,v 2 ) — d s (0, Vi, v 2 ). Thus, it is sufficient to assume that in the -y-direction fluctua- 
tions of p s are homogeneous. In this case, 

V(R)=e 2 J dvW(v)[d s (R,v)-d s (0,v)} (31) 

with W(v) = j dv + W(v + — v)W(v + + v ). The formula (31) together with Appendix B can 
be taken as a starting point to reproduce the results in LP00 using structure functions rather 
than spectra (compare with LP00). This is demonstrated in Appendix D. 

Substitution of the expression (26) into (29) gives 
P(R) = ^ j tfoiW>i) j dv 2 W(v 2 ) x 

X ^g~ a (PH+Pl2) _|_ g-a(P22+P2l) _ g-a(P21+Pl2) _ g-Ct(pll+P22)^ {32) 
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where we have used a shorthand notation p(Xi,Vj) = pij. Homogeneity in spatial directions 
dictates that after averaging the first term is equal to the second one and the third is equal 
to the fourth (e.g. for the density itself (pu+pi 2 ) = (P22 + P21) and (pi 2 +p 21 ) = {pn+Pu})- 

We shall now derive an approximate expression for X>(R), applicable for small separa- 
tions R. Let us rewrite eq. (32) as 

P(R) = ^ J dv 1 W(v 1 ) J dv 2 W(v 2 ) x 

X ^ e —a(ftu+Pi2) |~]_ _|_ g-a(P22+P2i-pn-pi2) _ e -a(p2i-pu) _ g -a(p22-pi2)j ^ (33) 

It is easy to see that the term in brackets depends only on the differences in density taken 
between two lines of sight at the same velocity, i.e. pu — p 2 \ and p 22 — p\ 2 - The point of 
rearranging the terms in this way is that for small separations R these differences are small, 
so we can retain only the leading order in the power series expansion for the corresponding 
exponentials 

P(R) ~ e 2 J d Vl W( Vl ) J dv 2 W(v 2 ) ( e -^ 1+ ^) [( Pll - p 21 ) (p 12 - p 22 )\) . (34) 

The structure of the resulting expression is as follows. The term in brackets is similar to 
the integrand in eq. (30) (indeed, ((p u - p 21 ) (p 12 - p 22 )) = d s (R,v 1 ,v 2 ) — d s (0,v 1 ,v 2 ) ) 
while the effect of the absorption is manifested in the exponential term which depends on 
the density distribution in velocity direction along the line of sight. As we see, the series 
expansion we have performed at small separation R is indeed an expansion that assumes 
a 2 [d s (R, v 1, v 2 ) — d s (0, v 1, 1*2)] to be small. Similar expansion is then applicable for any 
combination of structure functions that vanishes at R — and this opens a way of dealing 
with non-Gaussian statistics. 

Although there is a cross-correlation between the two terms in the integrand and sta- 
tistical averaging cannot, in a general case, be factorized, the correction for cross-correlation 
is of a higher order of smallness than the leading term. This can be seen by expanding the 
remaining exponential term in power series. We can illustrate this statement by specifying 
Gaussian statistical distribution for p s and performing explicit averaging which then gives 

V(R)~e 2 J dv 1 W(v 1 ) J dv 2 W(v 2 )e l i a2 ^ 1+pl2)2) [d s (R,v 1 ,v 2 ) - d s (0,v 1 ,v 2 ) + a 2 d 12 ] 

(35) 

The cross-correlation combination d 12 7 is vanishing as R — > and a 2 di 2 is small, as expected. 
For non-Gaussian statistics the irreducible higher order correlations appear. However it 



7 d 12 = [d s (R,v 1 ,v 2 ) - d s (0,v 1 ,v 2 ) + d s (R,v 1 ,v 1 )} [d s (R,v 1 ,v 2 ) - d s (0,v 1 ,v 2 ) + d s (R, v 2 , v 2 )] 
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is suggestive that these correlations provide vanishing contribution for R — > 0. Without 
providing a formal derivation for the general case, we propose that eq. (34) can be rewritten 

as 8 

T?(R)~e 2 j d Vl W( Vl ) J dv 2 W(v 2 )(e' a ^ + ^}[d s (K,v 1 ,v 2 )-d s (0,v l ,v 2 )] . (36) 
or when approximation of homogeneity in velocity direction is warranted 

P(R) ~ e 2 J d Vl W( Vl ) J dv 2 W(v 2 )(e- a ^ 1+pl2) ) [d s (R,v) - d s (0,v)} . (37) 

Expressions (36) and (37) show that for sufficiently small R the structure functions 
of intensity differ from the earlier studied case by the window function determined by the 
absorption W absorption = (e~ a(/>11+Pl2 )) . Two important conclusions directly follow from this 
observation. First, it is clear that if the window function determined by data slicing is much 
narrower than that given by the absorption we get results indistinguishable from the earlier 
studied case of no absorption. Second, in the case of integration over the whole line of sight, 
the results will be different from the earlier studied case because the window function is not 
equal to unity any more. The criterion when the absorption is not important for velocity 
studies is straightforward. Expression (37) transfers into (31) when W a b SO rption is of the order 
of unity over the range of scales studied. 



4. Scaling Regimes of the Emissivity Statistics 

In this section we shall study emissivity statistics in the regimes when it exhibits power- 
law scaling. Such behavior is possible at the small scales where the effect of the absorption 
is limited a 2 [d s (R, i>i, 1*2) — d s (0,vi,v 2 )] < 1 and eq. (37) provides a good approximation. 
We shall also restrict our consideration to the scales smaller than the cloud size R/S <C 1 
so that inhomogeneous boundary effects can be neglected. 



4.1. Integrated Lines 

The presence in eq. (37) of the window function determined by absorption, i.e. (e~ a( - Pll+p12 ^), 
changes substantially results compared to the case of no absorption discussed in LP00. In- 



8 A formal derivation of this in the case of Gaussian fluctuations is given in Appendix E. Our efforts 
above were directed to show that the structure of the expression does not depend on the assumption of the 
Gaussianity. 
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deed, the integration over the whole spectral line in the latter case removes the dependence 
of intensity statistics on velocity. The window function determined by absorption defines to 
what extend the integration over velocities is performed. 

Let us estimate the form of the window assuming Gaussian statistics and homogeneity 
of the density fluctuations 5p s = p s — (p s ). Then (see eq. (E3)) 

W absorpti on = e - a <P-( x i.^)> e - Q ^( x i.^)> e aa <^( x i.^)> e Q2 <^( x i.«a)> e -4*»( .«) . (38) 

For homogeneous statistics the variance of the fluctuations does not depend on velocity thus 
the corresponding factors are constant and affect only the normalization of the result. How- 
ever, their presence reveals limitations of the Gaussian approximation when the absorption 
is large. Indeed, a similar term appears in the mean line profile eq. (28) computed for Gaus- 
sian p s , namely, (J(X,Ui)) = ^ 1 - exp(-a(p s (X 1 , vj) + ^(<$p^(Xi, Ui))) . Clearly the 

answer is unphysical when a > 2(p s (Xi, v \)) / (8p^(X-i, v\j). The problem arises because the 
Gaussian fluctuations do not obey the constraint that the density p s is positive. Therefore 
for high absorption, negative density excursions, however rare, dominate the result. In view 
of this, we select in eq. (38) the factors which describe the variable part of the mean intensity 
profile, and write them in non-expanded form 

^absorption = (e~ aPs(Xl ) (e _aPs(Xl )e^~ ^ (Xl • t ' l)+<5p s ( x i ,V2)) e ~^-d s (0,v) 

This approximate formula does not suffer from the defect we have mentioned. The residual 
dependence on 5p s variance can be traced to the real effect of increase of intensity contrast 
between two points if absorption is present. However, this factor is constant in our treatment 
and we shall normalize it out. We summarize our considerations in the formula 

£>(R) oc J dv W(v) e-^ '") [d s (R, v) - d s {0, v)} , (40) 

where 

W(v) = J dv + W{v l )W{v 2 ){e~ ap ^' Vl) ){e~ ap ^' V2) ) . (41) 

The most important effect induced by absorption is an additional exponential down- weighting 
in the projection of the contribution from the points with large velocity separation v in a 
manner which itself depends on the turbulence statistics. For subsonic turbulence where the 
velocity variance is dominated by the thermal dispersion and is not scale dependent, the ef- 
fect of absorption in our approximation amounts to a constant and for statistics is equivalent 
to zero absorption and zero turbulent velocity case. However, one should still check whether 
the linearization of the kernel holds at the scales under study for a given absorption value. 
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The most important qualitative characteristic of the window is its width, which for 
absorption we shall define as velocity v ab at which 

a 2 d s {0,v ab ) = 1 . (42) 

In terms of the velocity-density decomposition of the PPV structure function, the product 
of the two windows arise ~ e -a 2 d v (o,v)/2 e -a 2 d p (o,v)/2 _ -q q ^ f ac tors act simultaneously but the 
one with the smallest width determines the gross effect. Asymptotic analysis of d s (0,v) in 
Appendix C gives for the absorption window widths 

v ab /D z (S) 1 ' 2 « (ap)^fc), m> 2/3 (l- 7 ) (43) 
Va b /D z (S) 1 ' 2 » (apy 1 , to < 2/3(1 -7) (44) 

where we have omitted numerical coefficients of order unity and estimated the mean density 
in PPV space as 

(p 3 )=pS/D,{S) 1/2 ■ (45) 
Absorption effects become negligible if v ab / D z (Sy/ 2 3> 1. 

Effect of Velocity Fluctuations 
Let us consider first the 5p s fluctuations that arise from velocity fluctuations only. Indeed, 
even if the underlying density is constant, random velocities do produce caustics that were 
shown in LPOO to be very important for the analysis. 

In Figure 1 we plot the correspondent kernel [d v (R,v) — d v (0,v)] as a function of v for 

2 - 

a fixed sample value of R and, simultaneously, the window function e~^ ds ^°' v ^ for the range 
of absorption amount characterized by ap s . The kernel is highly peaked at the zero velocity 
separation v — 0, has region of negative values, and approaches zero at large v. The window 
is unity for zero absorption but also peaks at v = as absorption increases. 

Integration of the velocity kernel over the whole line in the absence of absorption and 
W(v) = 1 gives zero, in concordance with the expectation of canceling all the velocity 
effects. In the presence of high absorption one might hope that down-weighting all nonzero 
v contribution will act as an effective thin slice reproducing thin slice power-law asymptotics 
of D(R), which is achieved when the width of the slice becomes small compared to rms 
velocity at the scale of the study R (LPOO). However we are restricted by the fact that 
our linearized approximation (34) is valid only when absorption is moderate a 2 d s (R, 0) < 1, 
while higher absorption (or larger scales at fixed a) will induce nonlinear modifications of 
the result. Thus we find thin slice scaling is never realized for to > 2/3. Instead, the new 
universal intermediate power-law scaling is exhibited for moderate absorption. Dependence 
of the 2D structure function on the underlying velocity index m is lost in this regime For 
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Fig. 1. — The velocity kernel d v (R,v) — d v (0,v) (thick solid line) and the density kernel 
d p (R,v) — d p (0,v) (dotted line) for R/S = 0.01 and m = 2/3 as a function of the line of 
sight velocity separation vj ' D l J 2 {S) normalized over velocity dispersion. The thin solid lines 
correspond to the window function exp(— a 2 d s (0, v)/2) and are drawn for the values a(p s ) 
from 0.3 to 30 (from top to bottom). For the displayed value of R, the linearized treatment 
of the kernels is applicable up to a{p s ) ~ v / 2(i?/S') (m ~ 2 ^ 4 6, the corresponding line is 
highlighted. 



m < 2/3 there is a range of scales when absorption can modify the scaling to that of the 
thin slice in which case sensitivity to m is present. 

Analytical estimates can provide some insight into the expected behavior of T> V (R) after 
integration as R — > 0. The value of the kernel at the peak, d v (R, 0) is given by the asymptotic 
expression (C4) (see also LP00 for asymptotics of the correlation functions in the velocity 
space) , 

d v (R,0) ~ 



V?r(m/4-l/2) 
r(m/4) 



pl-m/2 



(46) 
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where we are more interested in the R dependence of the expression than in the numerical 
prefactor. Incidently, this gives a criterion on the validity of the linear expansion of the 
velocity kernel. With dimensional quantities restored the condition 



sets the scale range where our treatment is applicable. 

Let us first consider the case m > 2/3. With the width of absorption induced window 
given by eq. (43), the condition (47) requires v 2 b > D Z (S)(R/ S) m , i.e velocity cutoff exceeds 
rms turbulent velocity at scale R and thus is insufficient to reproduce thin slice which require 
opposite relation to be true (see LPOO). 

This is, however, sufficient to suppress the negative large v tail in which case the domi- 
nant contribution to the integral scales as the height of the kernel peak times its width. The 
scaling of the width of the kernel peak with R can be readily obtained from the following 
consideration: If we integrate the kernel over v choosing fixed integration range dV but vary- 
ing R from large values to the small, then the projected result will scale with R according 
to thin slice asymptotics while D z (R) l l 2 > dV. The change to the thick slice behavior takes 
place when the condition is reversed. But formally in the thin slice regime the integral is 
determined by the zero-lag peak behavior only, while in the thick slice integration encom- 
passes most of the kernel. Therefore dV = D z (R) l l 2 , where D Z (R) is the z projection of the 
correlation tensor of velocity (see eq.(5), represents the required estimate of the width of the 
kernel. This scaling is valid both for the zero crossing point of the velocity kernel and the 
width of the density kernel that we consider later. 

The total scaling of the 2D emissivity structure function is thus T> V (R) ~ d v (R, 0)dV oc 
ffi-m/2 ^ ^ with dependence on m canceled out. The corresponding 2D spectra scales 



Numerical integration using eq. (35) confirms that indeed the velocity statistics exhibit 
a new intermediate asymptotic regime, which is characterized by a universal spectral index 
of —3. In Fig. 2 we show structure functions of intensity for two different indexes m = 2/3 
and m — 1. The range for the universal asymptotics is clearly seen. 

Fig. 2 also shows that for very small R the structure function grows faster than R, which 
entails spectra steeper than K~ 3 . Indeed, for sufficiently small R the exp(a 2 d v (0, v)) window 
does not cut out even the distant region where the kernel gets negative in which case we 
expect to see the thick slice asymptotics ~ R 1 + m / 2 (i. e x~ 3 ~ m l 2 for the spectrum, see LPOO) 
in accordance with our numerical calculations. Thermal effects (not included in numerical 
experiment) at scales where turbulence is subsonic D Z (R) < (3 lead to even faster fall-off of 




(47) 



as K~ 3 . 
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Fig. 2. — The structure function of intensity for m = 2/3 (right panel) and m — 1 (left 
panel) for a set of absorption coefficients a(p s ) = 0.2, 0.36, 0.7, 2.5 (from top to bottom). We 
use the general non-linear expressions from Appendix E, but plot the curves only up to R 
when linearized theory is still accurate. At the intermediate scales R the structure function 
tends to the universal regime ~ R which does not depend on the particular index m. The 
upper dashed line corresponds to the thick slice asymptotics for the velocity (spectral index 
—3 — m/2) and the lower dotted line correspond to the universal index —3 that does not 
depend on m. 

the structure function at small R and indeed make any coherent scaling of turbulent velocity 
hard to observe in emissivity statistics at that scales. 

For m < 2/3, on the other hand, absorption may lead to thin slice scaling oc R l - m / 2 
for some intermediate range of scales before R is so large as to cause nonlinear devia- 
tions. Indeed, in this case (ctp) v a b/D z (S) 1 ^ 2 and condition (47) can be rearranged to 
read D z (R)/v 2 b (R/ S) l ~ Zm l 2 < 1 which can be satisfied even for relatively high absorp- 
tion v 2 ah /D z (R) > 1 correspondent to thin slice velocity cutoff. For a given absorption 
parametrized by v a b we summarize the behaviour at the different scales in the Table 1. All 
four scaling regimes are present for m < 2/3 although how pronounced is actually thin slice 
scaling (4.1) depends on the extend of the scale range (v 2 b / D z (S)) l l m — (v 2 b / D z (S)) 2 / 2 ~ m , 
which grows with absorption strength if m < 2/3. This interval does not exist for m > 2/3. 
We note that the two-dimensional emissivity statistics of integrated lines T> V (R) restore m 
sensitivity for m < 2/3 where the line-of-sight statistics d v (0,v) is saturated at oc v 2 . The 
latter behavior is not surprising. Indeed, structure functions d v (0, v) cannot grow faster than 
v 2 and this explains the saturation. However, the information along the line of sight is not 
lost. The analysis of the spectra of fluctuations along the velocity coordinate performed in 
LP00 confirms this. 
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Scale Range 



Intensity Scaling 



Regime 



R/S<(P/D z (S)y/' 



velocity effects erased (subsonic regime) 

V(R) oc i^Wa (^icjfe s ii ce ) 

"D(-R) oc R 1 (intermediate scaling) 

V(R) oc R 1 -™/ 2 [thin slice) 



/ ,,2 \ 1/m 

« fe) 

/ „2 \ 1/m 
R I S < fe) 

(^)) Vm <^<l 

/ „2 \2/(2-m) 

(life) < «/S 




) 



2/(2-m) 



not a power law 



(strong absorption regime) 



Table 1: Scalings of structure functions of intensity fluctuations arising from velocity fluctu- 
ations for the power-law underlying 3D velocity statistics. In the strong absorption regime 
V does not follow a simple power-law. Thin slice regime does not exist for m > 2/3. D Z (S) 
is given by eq.(5). 

Effect of Density Fluctuations 
Density fluctuations modify the statistics of emissivity. Mathematically the density fluctua- 
tions are imprinted on the £ p part of the 3D correlation function. This correlation function 
depends both on the statistics of density and velocity. As we discussed earlier, density 
correlations may be short-wave or long-wave dominated. Within power-law model for over- 
density statistics they are described, respectively, by correlation £(r) = p 2 (r /r) 7 or structure 
d(r) = p 2 (ro/r) 7 functions, where 7 is positive in the former and negative in the latter case 
(see Appendix A). The amplitude of the density effect is encoded in the correlation scale r 
which has a meaning of the scale at which density inhomogeneities have an amplitude equal 
the mean density of the gas (i.e. bp ~ p for 7 > or p\ — P2 ~ p for 7 < 0). This gives rise to 
the amplitude factor {j-q/S) 1 in the asymptotic expressions of Appendix C (e.g. eq. (C5)). 

Contribution from density to PPV statistics combines linearly with the pure velocity 
term and an important question of which term dominates arises. Note, that both terms 
contribute to the absorption window at the same time in non-linear fashion. 

Table 2 summarizes at which scales the density inhomogeneity contribution dominates 
velocity effects in PPV space. For steep spectra, the density contribution dominates large 
scales, above characteristic scale determined by r , For shallow spectra it is dominant at 
small scales. Density inhomogeneities also dominate observed statistics at small scales where 
velocity is subsonic. 



What are the reasonable values for r ? We argue that the amplitude of density per- 
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Condition 



Case 7 < 



Case 7 > 



Eqns 



d p (R,0) > d v (R,0) : 



R>r 



R<r 



(C4, C5) 



d p (0,v) > 4(0, t;) : 
m > max [|, |(1 — 7)] 




v 2 <D g (S)(r /S) m (C8,C12) 



§ <m< |(l- 7 ) 
|(l-7) <m< § 



2 /3 not applicable (C8, C13) 



v 2 < D z (S)(r /S) (C9, C12) 

r /S>l (C9,C13) 



m < min [|, |(1 - 7)] 



Table 2: Range of the scales where the impact of density inhomogeneities to the PPV statistics 
exceeds the velocity contribution. For steep spectra 7 < 0, we need r < S for density effects 
to be seen. This is feasible in a model where high amplitude density inhomogeneities are 
saturated at scales r c < S. For the validity of the asymptotic scaling used, R <C S (or r c ) 
is required. When d p (0,v) > d v (0,v) the window function is determined by d p (0,v) (see 
eq. (42)). In addition, the window function is determined by density fluctuations for subsonic 
velocities. 

turbations of the scale of the cloud itself should not significantly exceed the mean density, 
otherwise the applicability of the notion of a cloud is suspect. For a shallow spectrum 7 > 
this means r < S. For a steep spectrum 7 < 0, r > S is likely, if a power-law distribution of 
fluctuation amplitudes extends all the way to the cloud scale S. In both cases (tq/S)"' < 1. 

An immediate conclusion from Table 2 is that for steep spectra 7 < density terms are 
expected to always be subdominant, since within a cloud R < S < r . However we must 
stress that this relies on an assumption that the turbulent density scaling extends all the 
way to the cloud size S, in which case high value of d(S) will indicate large inhomogeneities 
of the size of the whole cloud. But if the density structure function saturates at smaller 
separation r c < S, e.g. as d(r) = rf(oo)r _7 /(r -7 + r c 7 7 ), we may encounter a high amplitude 
d(oo) ^> p 2 without density at the scale of the cloud being significantly perturbed, since 
all inhomogeneities are restricted to smaller scales. As seen from the expansion d(r) ~ 
d(oo)(r/r c )~' y , r < r c , the effective correlation scale in this case is r = r c x [p 2 /ci(oo)] _1 / 7 < 
r c and in the range of scales r < R < r c one will observe the density effects. The physical 
picture, leading to such situation, is to have strong energy injection at scale smaller than a 
cloud size, with turbulent cascade establishing steep spectrum to smaller scales, while having 
energy sources only weakly correlated at larger scales. 
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The Table 2 shows that for a shallow spectrum the density term is important for PPV 
statistical descriptors in both R and v direction at small scales R < r and v 2 < D z (r ) (for 



The marked effect of absorption on intensity fluctuations in integrated lines is that the 
density contribution cannot always be recovered by increasing the slice thickness, as it is the 
case of no or small absorption. On the contrary, we are able to observe the density impact on 
the emissivity only in restricted scale range. The reason for this is that for a wide range of 
scales above the thermal dispersion scale the velocity term is not suppressed in the presence 
of the absorption even when we integrate over the whole line. 

We shall now discuss the how the density contribution scales, depending on separation 
R and conditions on 7 and m. Whether density term will dominate the overall PPV statistics 
depends on the amplitude of density inhomogeneities given by r according to Table 2. Our 
results are summarized below in Table 3. 

Scaling estimations for density are obtained in a similar way to the case of pure veloc- 
ity term in the previous section. The main ingredients are density kernel and absorption 
windows, exemplified in Figure 1. The complication is that the absorption window is de- 
termined by d s (0,v), with contributions from both velocity and density giving rise to two 
multiplicative windows (see equation (42)). The one with the smallest width of the two is 
the most important one. We shall explicitly note when the details of the absorption window 
are important. 

The peak of the density kernel now scales as R}-i- m / 2 ( se e equation C5) while its width 
is dV ~ D^R) 1 / 2 , as we discussed before. Integrating over the whole line in the absence 
of absorption gives an estimate T>(R) oc R 1 -~i- m l 2 R m l 2 ~ R}-i (this refers to a correlation, 
rather than to a structure function of emissivity if 7 > 1). which is exactly the thick slice 
regime of LPOO. If absorption is not too strong, v a b > D^R) 1 ^ 2 , it will not affect this result. 
This is different from pure velocity term, since the velocity term exactly vanishes when 
integrated over the whole velocity range in the absence of absorption, and small absorption 
leads to qualitatively new intermediate "universal" scaling. 

The condition of weak absorption in the presence of density, that replaces the one in 
(47), is given by 



Using (43,44) to express absorption via v a b, we find familiar constraints when absorption 
effects are not strongly non-linear 



m > 2/3). 




(48) 



v 2 ab > D z (S)(R/Sr~D z (R), 



m > 2/3 (I-7) 
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Scale Range 



Intensity Scaling Regime 



R/S < {P/D^S)) 1 /™ V(R) oc R 1 -"/ (subsonic) 

/ 2 \ l/ m 

R/S < (^y) oc i? 1 -^ (thick slice) 

(^)) 1/m <^<(z&)^ W) « ^- /2 (to slice) 

/ „2 \2/(2-m) 

( d^5) ) < -R/'S' n °t a power law (strong absorption) 

Table 3: Scalings of structure functions of intensity fluctuations arising from density inho- 
mogeneities for the power-law underlying 3D velocity statistics. Thin slice behaviour require 
7 < 1 — 3m/ 2. In the strong absorption regime T> does not follow a simple power-law. 

v 2 ab > D Z (S)(R/Sy^~ m / 2 , m< 2/3 (l- 7 ) (50) 



The first case, m > 2/3(1 — 7) is applicable for the most interesting values of m if the 
density spectrum is shallow, 7 > 0. Here, linearity and thick slice conditions coincide and 
one cannot achieve thin slice regime before the absorption effects become nonlinear. 

The second case, m < 2/3(1 — 7) has wide range of applicability for steep spectra 7 < 0. 
Here linearity condition is less stringent and in the range of scales (vl b /D z (r c )) 1 / m < R/r c < 
(vl b / ' D z (r c )) l l l ~' 1 ~ rn l 2 absorption effects modify the emissivity scaling towards the thin slope 
jj>i- 7 -m/2^ w j 1 ji e nonlinear effects come into play only at larger scales. 9 . 

To summarize: For shallow spectra one expects line-integrated statistics for density 
perturbations to follow thick slice scaling until absorptions effects become too strong. Density 
terms are visible in the overall emissivity statistics at R < r . For steep spectra, density 
contribution may exhibit, in addition to thick slice behaviour at small scales, the thin slice 
scaling at intermediate scales. However, density contribution can be dominant in emissivity 
statistics only when R > tq, which require high amplitude density inhomogeneities saturated 
at scales smaller than a cloud size. 

Table 3 and Fig. 3 illustrate how the statistics of emissivity changes with the change of 
the spectral index of density and velocity. 

Can we ever observe density fluctuations in the absorbing gas outside the ranges set by 



9 To be consistent, velocity effects of Table 1 have to be reevaluated in this case due to larger contribution 
from density perturbations to the absorption. 
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Fig. 3. — Left panel: The structure function of intensity for Kolmogorov density and Kol- 
mogorov velocity. Right panel: The structure function of intensity for shallow density 
7 = 1/2 Kolmogorov velocity m = 2/3. Absorption parameter for the lines shown varies 
in the range a(p s )(r / S) 1 ^ 2 = 0.1 (top) to 1 (bottom). Scaling solutions for both cases are 
given in heavier line and labeled. The calculations are done using unexpanded formula given 
by eq. (E4). The non-linear effects get important as R increases. 

correlation scales, for example when turbulent velocities are small ? While one decreases the 
amplitude of supersonic turbulence the absorption effects increase due to atoms concentrating 
at similar velocities which is expressed by PPV density p s pS/ y/D z (S) + 2(3 growth, 
making recovering turbulence scaling more difficult. However velocity effects disappear as 
turbulence becomes subsonic at all scales D Z (S) <C j3. The power law index of density 
inhomogeneities can then be obtained over the range of scales if apS/ \ffi <C 1. Under similar 
conditions we can see density fluctuations at small scales r < /thermal > where the turbulence 
is subsonic even if it is supersonic at large scales. The scale r thcrmal can be determined as 

^(rthermai) = fl- 
it may look somewhat counterintuitive that even in the absence of velocity we never see 
a thin skin of the cloud exhibiting 2D spatial 10 slice of density fluctuations. As absorption is 
proportional to density in our model, the optical depth changes with the change of density. 
As the result the emission fluctuations are not dominated by the atoms within the thin skin 
depth from the cloud surface. The origin of the effect can be traced to original equations 
of radiative transfer (see eq. (26)). For instance, in the absence of velocity effects p s is 
proportional to the integral over the whole line of sight (see eq. (25)). 



10 We may note that a thin velocity slice is very different from a thin slice in space. Integration over the 
entire spatial extend of the cloud is performed while the velocity slice is defined. 
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4.2. Thin Velocity Slices of Spectral Line Data 

Within the VCA technique introduced in LPOO the change of the spectral index of 
emission fluctuations with the thickness of the channel maps was employed to separate the 
contributions from density and velocity. In particular, to see the thin slice scaling that 
reflects the statistics of turbulent velocities, the effective window of the velocity slice should 
be narrower than the turbulent velocity dispersion at the scale of the study. 

The thermal broadening (i.e. /3 factor in the previous section) of the line can be con- 
sidered as a convolution with the window W s i ice describing the channel slicing adopted in 
the data analysis. The LPOO study showed that adopting the slice thickness less than the 
thermal width 8V < (3 does not provide new information about turbulence 11 . In our discus- 
sion we shall treat thermal effects as adding (in quadrature) the minimum width to the data 
slice, symbolically 5V 2 + (3 

Absorption results in appearance of the second window function W a b SO rption- Two cases 
can take place. 

If v 2 b < 5V 2 + (3 the regime is not different from that discussed in the earlier subsection. 
Indeed, it is the absorption that limits the integration and the results of integration over the 
slice thickness and over the whole emission line should be similar. According to previous 
discussion, absorption smoothing reveals the underlying velocity statistics only for m < 2/3. 
Since the thermal broadening is given by the properties of the media under study and 
therefore is fixed, we get the following picture. If absorption is strong v 2 b < (3 making slices 
thinner will not provide any additional information about turbulence. This is the case when 
the velocity information recovered using the VCA is limited if any. 

In case v 2 b > (3 one can get information about turbulence by making the slice thinner, 
i.e. decreasing SV to 5V 2 < v 2 b — (3. Choosing the thickness of the velocity slice to ensure 
that v 2 b > 5V 2 + (3 we formally return to the regime very similar to that considered in LPOO. 
Indeed, if the window function defined by the absorption is unimportant we return to the 
case mathematically similar to that in LPOO. Naturally our results obtained in this regime 
are identical to those in LPOO. Fig. 4 shows what is happening when we choose thin slicing 
for self absorbing data. It is easy to see that in doing so we get back to the regime of thin 
slices discussed in LPOO. 

As the result we can extend the VCA for self-absorbing clouds and provide observers 



n It docs provide the information about the gas temperature, however. This point will be elaborated 
elsewhere. 
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with a following recipe: 

1. If the integrated line emissivity shows a universal spectrum with n = —3, this testify 
that the density statistics is steep. In this case thin slice regime would provide one with a 
spectrum, which power-law index is —3 + m/2. 

2. If the integrated line emissivity shows a spectrum which is shallower than —3, this index 
should be identified with the power-law index n of the underlying 3D density field. The thin 
slice regime in this case provides the index n + m/2. 

It is clear that in the case of shallow density the VCA recovers both 3D velocity and 
density spectra, while in the case of steep density the VCA would normally recover just 
velocity spectral index. Potentially, one may attempt to remove the universal statistics in 
this case, but this may be very challenging in practical terms. It is also evident that for 
subsonic turbulence thin slice regime is not attainable. 

5. Accomplishments and Qualitative Understanding of Results 

Let us start with an explanation what has been done so far in different parts of the 
present paper. We started with introducing radiative transfer in §3 and obtaining general 
expressions for the statistics of intensity in the presence of the radiative transfer. We showed 
that absorption plays the role of an additional window function in those expressions. Then in 
§4 we discuss the two limiting and observationally important cases. One is related to the way 
how most present day observations of fluctuations are performed, namely, through measuring 
fluctuations of the integrated intensity. In this situation we show that the information about 
the underlying statistics in most cases is lost. We show, however, that the information can 
be recovered if fluctuations in velocity slices are analyzed. These are two major results of 
the paper. 

A good deal of the paper and most of the Appendixes deal with the statistics of the 
fluctuations without absorption. This was necessary because the earlier machinery was 
formulated in terms of power spectra, while the transformation from the PPV power spectra 
to correlation functions happen to be not trivial. In terms of our earlier work, this not only 
provides an independent check of our results in LPOO, but provides a very important extension 
of the technique. For instance, we obtain the predictions for the statistics of fluctuations 
along the velocity direction. These fluctuations can be studied using instruments with good 
spectral resolution even if the spatial resolution of the instruments is poor. 

With the mathematical machinery in the paper being sufficiently involved, the issue 
of qualitative understanding of our results is of a particular interest. Especially this is 
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important for our results that may look somewhat counterintuitive. 

Consider first the formulae for thin slices. The spectrum of intensity in a thin slice 
gets shallower as the underlying velocity get steeper. To understand this effect consider 
turbulence in incompressible optically thin media. The intensity of emission in a slice of 
data is proportional to the number of atoms per the velocity interval given by the thickness 
of the data slice. Thin slice means that the velocity dispersion at the scale of study is larger 
than the thickness of a slice. The increase of the velocity dispersion at a particular scales 
means that less and less energy is being emitted within the velocity interval that defines 
the slice. As the result the image of the eddy gets fainter. In other words, the larger is the 
dispersion at the scale of the study the less intensity is registered at this scale within the 
thin slice of spectral data. This means that steep velocity spectra that correspond to the 
flow with more energy at large scales should produce intensity distribution within thin slice 
for which the more brightness will be at small scales. This is exactly what our formulae 
predict for thin slices (see also LPOO). 

The result above gets obvious when one recalls that the largest intensities within thin 
slices are expected from the regions that are the least perturbed by velocities. If density 
variations are also present they modify the result. When the amplitude of density pertur- 
bation becomes larger than the mean density, both the statistics of density and velocity are 
imprinted in thin slices. For small scale asymptotics that we are mostly interested in the 
paper, this happens, however, only when the density spectrum is shallow, i.e. dominated by 
fluctuations at small scales. 

Understanding of the results for integrated spectral line in the presence of absorption 
is a bit more involved. Again, it is advisable to consider incompressible flows first. If no 
absorption present the integrated spectral line images are not affected by velocity fluctu- 
ations. Indeed, the velocity field just redistributes atoms along the line of sight and this 
redistribution cannot affect the total intensity in the absence of absorption. 

If absorption is present, the fact that velocity redistributes atoms along the line of sight 
causes fluctuations of the integrated intensity. Two effects, however, are present simulta- 
neously as the velocity field at a particular scale spreads atoms in PPV space. It is easy 
to see that the most absorption is expected is when the atoms have the same velocities. 
Therefore velocity dispersion provided by turbulence decrease the absorption and therefore 
increases the signal. However, the decrease of the optical thickness makes the media to be- 
have more like its optically thin counterpart, in which the velocity fluctuations are averaged 
out. In other words, while the mean level of intensity increases as the result of turbulence, 
the contrast of the fluctuations caused by the velocity fluctuations decreases. The range of 
the universal asymptotics K~ 3 is the range at which the two effects exactly compensate each 
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other. 

Density fluctuations get important only when their amplitude gets larger than the mean 
density. This is quite analogous to the case of thin slices discussed above. The outcome is 
also analogous, namely, for the small scale asymptotics only density with shallow spectrum 
is important. Density with steep spectrum may affect the large scale asymptotics. 

6. Discussion 

6.1. Description of turbulence in absorbing media 

In LPOO we provided the statistical description of spectral line data cubes obtained via 
observations of turbulent media. In the present paper we generalized this description for the 
presence of the absorption effects. 

If we compare our present treatment of PPV statistics with that in LPOO we see that 
it differs in terms of the statistical tools used. As the absorption takes place in real and not 
Fourier space, the use of structure functions is advantageous. In Appendix D we re-derive 
LPOO results for optically thin media using structure functions. Our present treatment 
extends the VCA technique for studies of turbulent absorbing media. How good is our 
treatment of absorption? 

Above we used ID model of radiative transfer. We would, however, argue that it is 
adequate for our purposes. Indeed, our aim is to study velocity and density fluctuations 
before the effects of optical depth distort the underlying power-law spectra. On the contrary, 
when the 3D radiative transfer gets important we expect the non-linear flattening of the 
observed emissivity spectra. 

For the sake of simplicity in our treatment we made major calculations assuming that the 
fluctuations within PPV data cubes are Gaussian. However, we showed that the structure of 
our expressions does not depend on this assumption. Whatever is the statistics of emissivity 
fluctuations in PPV space the absorption introduces a new window function that truncates 
the velocity range over which the integration is provided. Possible changes in the window 
function arising from non-Gaussian fluctuations are expected to result in the unimportant 
changes of the numerical prefactors of the expressions. 

We can test our result in the limiting cases of high and low absorption. In the latter 
case our formulae can be expanded and we get the LPOO result (see also Appendix C) for 
a — > 0. For high absorption, it is only very small fluctuation that are still in the linear 
regime. Those small fluctuations are likely to be Gaussian. 
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There are a number of issues that are relevant to studies of turbulence in both absorbing 
and not absorbing media. One of them is importance of velocity-density correlations. While 
the issue, in general, require further research, the case of velocity-density correlations has 
been studied analytically in LP00 and numerically (Lazarian et al. 2001, Esquivel et al. 
2003). Our discussion above reveals that when the effective thickness of the slice is less than 
window function given by absorption W 'absorption the effects of absorption are unimportant. 
For such slices, if they are thin, the effect of velocity-density correlations was shown to be 
marginal. 

For the sake of simplicity both in our present paper and in LP00 we assumed that the 
emissivity is proportional to the first power of density. If the dependence is non-linear, e.g. 
proportional to n 2 as in the case of H a emission, for small amplitude perturbations it is 
possible to show in the linear terms dominate and therefore our results above are applicable 
(see discussion in Cho & Lazarian 2002b) 12 . 

An additional simplification shared by this work and LP00 is that the thermal velocity of 
media is assumed constant. As discussed above thermal velocities affect the effective window 
function over which the turbulent velocities are smeared. As the result if observations sample 
regions with different temperatures of gas, the thermal smearing will provide higher weighting 
to the regions where thermal velocity is less than the slice thickness. This stays true in the 
presence of absorption. 

Another non-linearity disregarded in the treatment above as well as in LP00 is related to 
the non-linearities in the Galactic rotation curve. The effects of the rotation curve shear were 
studied in Esquivel et al (2003). There it was shown that effects of the shear can be ignored 
unless the galactic shear and the shear provided by the eddies are comparable. The latter 
is a very artificial situation as the large scale shear in most cases produces turbulence. For 
Kolmogorov turbulence the local shear vjl scales as Z~ 2 / 3 . Magneto-rotational instability 
(see Chandrasekhar 1967, Balbus & Hawley 1995) ensures that the global galactic shear 
creates turbulence. If the shear provides marginal influence, the higher order corrections 
provided by the non-linearity of the shear are also negligible. 

Our treatment above assumes that the absorption happens within the emitting gas. 
What will happen when the radiation source is outside the turbulent cloud? If turbulence is 



12 One may wonder whether complications arise at large scales. We believe that this is unlikely for steep 
densities. For instance, for H Q the emissivity is proportional to density squared. The correlation functions 
of density 2((5pi) 2 (5p2) 2 } correspond to a steeper spectra and therefore their contribution within thin slices 
is sub-dominant. 
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studied using absorption lines the registered intensity is 



I = Iq exp 



a(z)n(z)<f>(z)dz 



(51) 



where Jo is the intensity coming from the external source. The easiest way is to correlate 
logarithms of I /Iq. For constant a these quantities are evidently proportional to the integrals 
of emissivities and the analysis in LPOO is directly applicable. 

Is it always true that we have to deal with velocity fluctuations? If the velocity turbu- 
lence gets damped, while density statistical structure persists, e.g. in the form of entropy 
fluctuations, the effects of velocity modification get marginal 13 . 



6.2. Interpretation of Observations 

In a number of instances the power-law spectra of interstellar turbulence were reported 
for the spectral line data. They include HI (Green 1993, Stanimirovic et al. 1999) CO 
in molecular clouds (Stutzki 1998), H a in Reynolds layer (Minter & Spangler 1996), HI 
absorption (Deshpande, Dwarakanath & Goss 2000). In all the cases the reported index was 
close to —2.7, which is really amazing as in all those cases very different quantities were 
measured 14 . In comparison, fluctuations of synchrotron emission and starlight polarization 
show very different indexes, although underlying turbulence is likely to be Kolmogorov (Cho 
& Lazarian 2002). 

Results of the earlier applications of VCA to the HI data can be summarized as follows: 
I), the observed scaling of 2D power spectra within thin slices is consistent with arising 
from velocity caustics; II). the density power spectrum is steep at least in the case of Small 
Magellanic Cloud data; III). HI data is consistent with the velocity spectral index being 
approximately Kolmogorov. 



13 An important case of this regime is the viscosity damped MHD turbulence reported in Cho, Lazarian 
& Vishniac (2002c). This new type of MHD turbulence produces small scale magnetic fluctuations through 
stretching of magnetic field by marginally damped eddies ( Lazarian, Vishniac & Cho 2003). The slowly 
evolving fluctuations of magnetic field result in the fluctuations of density that are essentially stationary 
(Cho & Lazarian 2003ab). 

14 The spectral index obtained in Stanimirovic et al. (1999) was somewhat steeper, namely, —3. However, 
it was shown in Lazarian & Stanimirovic (2001) that this difference stemmed from the difference in the slice 
thickness. As slice thickness was reduced the index became w — 2.8 over the whole range of scales under 
study. 
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The interpretation of H a data may be done along the same line of reasoning. Earlier we 
argued that non-linearities in the H Q emissivities are not important for studies of small-scale 
fluctuations. This would also indicate the correspondence with Kolmogorov spectrum of 
turbulence of velocity. 

Consider now HI absorption. Instead of the over-density the absorption measures deal 
with the ratio of density over temperature n(x)/T(x). Both quantities fluctuate and a priori 
one expects the spectrum to differ from that of HI emission. However, Deshpande et al. 
(2000) found that the index of the absorption spectrum in slices of their data is similar 
to that Green (1993), but over the range from 3 pc to 0.07 pc. However, at scales of 
approximately 0.2 pc a change of the index of the power spectrum was observed. Therefore 
the observed power law covers approximately one order of magnitude scales from 3 pc to 
0.2 pc. If at those scales turbulent velocities are larger than thermal velocities, the slices 
used may be thin and the coincidence of the spectral index in the emission and absorption 
would mean the continuation of the turbulent cascade to small scales. This, however, may 
be a bit surprising as the velocity turbulence at those scales is likely to be damped through 
viscosity by neutrals unless the turbulence is generated at very small scales. 

To determine what is going on Deshpande et al. (2000) used VCA prescriptions to 
evaluate the effects of the velocity and did not detect the change of the spectral index as he 
changed the thickness of the data slice. Deshpande et al. (2000) could not see the changes 
of the spectral index and this can be interpreted that the fluctuations arise exclusively from 
density variations with the shallow spectral index (~ —2.7). 

If further research confirms that the density spectrum is indeed shallow this may be 
the signature of the magnetic energy cascade simulated in Cho, Lazarian & Vishniac (2002), 
Cho & Lazarian (2003) and described in Lazarian, Vishniac & Cho (2003). In this new 
regime of turbulence fluctuations of density arise from fluctuations of magnetic field, while 
velocity fluctuations are marginal. The expected power spectrum of density scales k~ 3 and 
is roughly consistent with the observations. Note, that such a shallow spectrum of density 
according to Deshpande (1999) can explain the detected mysterious tiny scale structures in 
HI (see Heiles 1998). It is obvious that a further study of the statistics of HI fluctuations at 
small scales is necessary. Synergy of different techniques can be called for. For instance, 
modified velocity centroids (Lazarian & Esquivel 2003) can provide an independent test 
whether velocity fluctuations are important. 

All the data above was obtained by studies of fluctuations within velocity slices of 
spectroscopic data. On the contrary, CO data by Stutzki et al.(1998) present the power- 
spectrum of fluctuations obtained via integrating the intensity over the whole emission line. 
According to Falgarone et al. (1998) the turbulent velocities even at the smallest scales are 
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still supersonic. Thus we expect to see a velocity modification of the spectrum. Stutzki 
et al. (f998) get the index of order —2.8. Accepting possible criticism that in the case of 
HI we assumed that the deviation of —2.7 from —3 is meaningful and informative, we may 
speculate that the measured CO spectral index are close to the universal asymptotics of 
—3 predicted in the paper above. We may speculate that the flattening of the spectrum 
compared to the expected —3 may be caused in the case of CO by instrumental noise (this 
possibility was mentioned to us by our referee). 

Alternatively, one expects to see the power index n + m/2 if the density spectral index 
n > —3. If we want to account for —2.8 index we would have to accept not only that the 
density is shallow, but also that the velocity is more shallow that the Kolmogorov value. 
More studies utilising our present theoretical insight are clearly necessary. 

What are the goals that can be achieved by studying properties of interstellar turbu- 
lence? Establishing turbulence spectra by applying VCA to different emission and absorption 
lines can help to answer fundamental questions of the interstellar physics. For instance, it 
can answer the question whether the turbulence in molecular clouds is a part of a bigger cas- 
cade from the large scales or it is locally generated. It can also answer the question to what 
extend the self gravity is important in determining the spectra of fluctuations at different 
scales. 

In addition, VCA states that the intensity of fluctuations measured in a velocity chan- 
nel depends on whether the channel width is larger or smaller than the thermal velocity 
dispersion of the atoms. This potentially allows to find the distribution of temperatures of 
emitters along the line of sight. This may provide an independent testing of theories of the 
ISM thermal structure and answer the question of the role of thermal instability for the ISM 
(see Heiles & Troland 2003). 

6.3. Range of Applicability and Numerical Testing 

The VCA technique is the most efficient for restoring the underlying spectral indexes 
when the turbulence obeys the power-law and the turbulent velocities are larger than the 
thermal velocities of atoms. While the first requirement is usually fulfilled in interstellar 
turbulence (see Stanimirovic & Lazarian 2001), numerical simulations provide a limited 
inertial range if any. This presents problems when the testing is involved. For instance, in 
Lazarian et al. (2001) and Esquivel et al. (2003) the spectral slope was corrected to fit a 
power-law. While the power-law requirement is not an essential constraint for most of the 
interstellar situations, the ratio of the turbulent to thermal velocity may be a constraint. To 
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use the VCA for the subsonic turbulence one should use the heavier species, which thermal 
velocities are less than that of turbulence if the thin slice asymptotics is intended to be 
observed. Note, that in case of turbulence that is supersonic at large scales this potentially 
limits the range of the small scales which can be studied through the VCA technique. Indeed, 
at some small scale the turbulence gets subsonic and the slope of the emissivity fluctuations 
within a slice should change 15 . 

Another problem of numerical testing is the appearance of shot noise related to the 
limited number of data points available through simulations. Unless understood, this noise 
produces a bias in the power spectral indexes obtained through synthetic observations based 
on the numerically generated data cubes. This problem was analyzed in Esquivel et al. 
(2003), where it was concluded that it is not important for actual observations where the 
number of emitters is essentially infinite. When absorption is taken into account our prelim- 
inary analysis shows that the testing is becoming even more challenging and requires piling 
up many numerical data cubes in order to obtain adequate statistics. 

The non-trivial issues above may cause confusion. For instance, in Miville-Deschenes et 
al. (2003a) the utility of the VCA was questioned. We have to reiterate that the effects that 
the authors encountered, e.g. shot noise, are related to the synthetic data they used and 
not to the actual observational data. None of the problems they discuss was present when 
supersonic turbulence of SMC was studied in Stanimirovic & Lazarian (2001). 

In another paper Miville-Deschenes et al. (2003b) mention limitations of the VCA 
related to the thermal motions of atoms. It is pointed out there that the thermal line width 
of atoms can entail asymptotic behavior that can be mistakenly identified with the actual 
saturation of the power-law index in a thin slice. In fact, the observed saturation in the two 
cases is different. If the thermal line width gets larger than the thickness of the slice, further 
decrease in the slice thickness does not change the shape of the slice emissivity spectrum. It 
consists of two parts, a shallow part corresponding to the scales at which turbulent velocity 
is larger than v t hermai (thin slice) and steeper part corresponding to the scales at which the 
turbulent velocity is larger than v t hermai (thick slice). On the contrary, if the turbulence is 
supersonic on the scales under the study, the change of the slice thickness makes the change 
from thin to thick slice to move to smaller and smaller scales. In other words, the slice is thin 
at large scales and thick at small scales with the border between these two regime moving 
as the thickness of the slice changes, provided that the experimentally determined thickness 



15 This change can be used to determine at what scale turbulent and thermal velocity are equal. If the 
dispersion of turbulent velocities is known at the large scales, this allows to measure the temperature of 
gas. If, alternatively, the temperature of gas is known, this allows to calibrate the power spectra in terms of 
absolute values of turbulent velocities. 
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of the slice is larger than the thermal broadening of the line 16 . 

Such a study has not been done in Miville-Deschenes et al. (2003b). Neither they give 
estimates of the thermal and turbulent velocities of gas in the region Ursa Major Galactic 
cirrus they studied. Instead they pointed out to the discrepancy between the spectral index 
obtained using the velocity centroid analysis and those the VCA would give if the asymptotic 
regime were not due to thermal broadening. By itself, this shows important synergy of 
different techniques. Unfortunately, the authors did not check whether their centroids are 
dominated by velocity or density. The application of the criterion obtained in LE03 to the 
data seems essential for answering this question. High degree of correlation between the 
averaged density and the averaged centroids maps as well as a close similarity of the spectra 
deduced from averaged density and the centroid maps calls for further testings of what the 
velocity centroids measure velocity in the particular case (see LE03). A study in Esquivel 
& Lazarian (2004) testifies that velocity centroids are dominated by density if magnetic 
turbulence has Mach numbers higher that 2. This is also suggestive that centroids obtained 
for the hypersonic turbulence in HI are dominated by density. 

We would like to stress that the importance of quantitative understanding of the statis- 
tics of spectral line data cubes goes far beyond the particular recipes of the VCA technique. 
LP00 and this paper provides for the first time the machinery suitable for dealing with the 
statistics of spectral line data. This machinery should provide advances not only when anal- 
ysis of spectra in channels is involved, but also suggests new techniques. For instance, our 
treatment of velocity correlations in v direction allows a new way of measuring the spectral 
index of velocity fluctuations. 

7. Summary 

In this paper we have studied effects of absorption on the statistics of spectral line data. 
In particular, we discuss how the 3D velocity and density statistics of turbulence can be 
recovered from the spectral line data. We subjected to scrutiny two regimes of studying 
turbulence, namely, (a) when integrated spectral line intensity is used and (b) when the 
channel maps are used. For our treatment we assumed the power-law underlying velocity 
and density statistics. Our results are summarized below. 

I. For integrated spectral line data we obtained: 



16 The good news is that by varying the thickness of the velocity slices on can find the ratio of the thermal 
and turbulent velocities (LP00). 



-36- 



1. Absorption introduces non-linear effects that distort the power law behavior of mea- 
sured intensity. The non-linear effects are marginal when fluctuations are studied at suffi- 
ciently small scales. 

2. While in the absence of absorption the velocity effects were washed out for the 
integrated spectral data, absorption makes velocity modification prominent. The recovery of 
the underlying velocity spectrum from integrated spectral data is non-trivial and sometimes 
impossible. For a range of scales the resulting power spectrum of intensity is caused by 
velocity fluctuations but does not depend on the spectral index of velocity fluctuations. 

3. Density spectra can be recovered in the limit of small scale asymptotics that we mostly 
deal here only when the density has a shallow spectrum, i.e. is dominated by fluctuations 
at small scales, or at scales where turbulence gets subsonic. 

II. For spectral data in channel maps we obtained: 

1. If slices are sufficiently thin we recover thin slice asymptotics obtained earlier in 
LPOO. However, the choice of slice thickness is limited by absorption, which introduces the 
minimal thickness of the slice for which non-linearity is negligible. 

III. We incorporated our results into the VCA technique intended for the observational 
studies of turbulence. The detailed description of the spectral line data cube statistics that 
we obtained should enrich the set of tools available for advancing quantitative insight into 
interstellar turbulence. 

Part of this work was initiated while one of us (AL) enjoyed the stimulating environment 
of the Astrophysical Turbulence Program 2000 at the Institute for Theoretical Physics at 
the University of California at Santa Barbara. AL also acknowledges the Visiting Professor 
position at Ecole Normale that enabled some of the discussions helpful for this research. 
Research of AL is supported by the NSF Grant AST-0307869 and by NSF Center for Mag- 
netic Self-Organization in Laboratory and Astrophysical Plasmas, while research of DP is 
supported by the Natural Sciences and Engineering Research Council of Canada. 

A. Statistics of velocity and density in XYZ coordinates 

Statistics of a finite-size cloud. If we deal with a cloud of size S the turbulence has 
a maximal scale. In the absence of additional physics the size S also serves as structure 
function cutoff scale D(r) ~ D(oo)r m /(r m + S m ). We notice that if this law is adopted for 
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DiL^r), z-component structure function in the case of solenoidal turbulence is 
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We shall further assume that the turbulent velocity field is Gaussian, in which case its 
properties are fully determined by the two point probability distribution function P(ui,u 2 ), 
where U\ = ■u 2 (x 1 ), u 2 — w z (x 2 ). Using variables u — U\ — u 2 , u + = (ui + u 2 )/2 this function 
is conveniently expressed (we write down the distribution of z-component of the velocity field 
■u(x) only) as 
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Validity of eq. (A3) is immediately checked by observation that u, u + are uncorrelated 
(uu + ) ~ (uf) — = Gaussian quantities with dispersions 
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and the Jacobian of transformation between u, u + and Ui,u 2 is unity. 

Shallow density spectrum 
u i — M z( x i) ; w 2 = m 2 (x 2 ) When statistical properties of the density fluctuations dominated 
by short wavelengths, we use power-law correlation functions of over- density: 



*(r) = (p) 2 (l + 



r \T 



7 > . 



(A4) 



Note, that the power-law part of the correlation corresponds to the 3D power-law spectrum 
oc k n , where n = 7 — 3 which for 7 > means that n is shallow, i.e. less than —3. 

Steep density spectrum 

To describe the density statistics for steep (i.e. n < —3) power spectrum one should use 
structure function description. For 7 < we write 



e(r) = -d(oo)[l-d(r)Moo) 



(A5) 



where d(r) is a structure function of density given by Eq. (3) 
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Real world structure functions do not grow infinitely and therefore we have to introduce 
a cutoff at some large scale. If the cut-off happens at r c then for the long-wave dominated 
turbulence (7 < 0) 



r 



r -7 + r c 
and the correlation function is 

2 r~T + r. 
For sufficiently small r Cr c eq. (A7) gives 



d(r) = d(oo) - ■ _^ , (A6) 



«0 = ^7=^=r • ( A7 ) 



£(r)«±d(<x))(l-[r/r c p) (A8) 

Note that in eq. (A8) r c [(p) 2 /d(oo)]~^ plays the role of r in eq. (A4). 

When using both correlation and structure functions one essential difference should be 
kept in mind: while the value of the structure function at some scale r*, d(r*), is determined 
by the power of fluctuations at smaller scales r < r*, the value of the correlation function £(r*) 
reflects the integral power at scales r > r*. For the power-law statistics, the correspondent 
integral contributions are dominated by the power at scale r* itself, localizing the information 
provided by £(r*) and d(r*) values. 



B. Statistics of density in PPV 



In this Appendix we calculate some statistical properties of the velocity space density 

(Bl) 



p s (X,v) = f dz p(x)0 w (x), 
Jo 



where 0„(x) is given by eq. (13) 

and random fields p(x) and -u(x) over distribution of which we average are assumed to 
be uncorrelated. We assume turbulent velocity field to be described by the Gaussian two- 
point probability distribution function given in eq. (A3). We shall also consider statistical 
properties of the density distribution in the galactic coordinates to be homogeneous: 

The mean PPV density is given by 



(p s (Xi,ui)) 



dzi (p(xi)) (0„i(xi)) 



(tt [D z (oc) + 2P]) 1/2 Jo 



dz\ exp 



( v l — v gal,l) 2 

D z (oo) + 2(3 



;b2) 
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and the two point correlation function is 



(p s (X 1 ,v 1 )p s (X 2 ,v 2 )) = dzi dz 2 (p(xi)p(x 2 )) (0„i(xi)0„ 2 (x 2 )) 

Jo Jo 



2tt /_ 



= — dz 



X 



S-\z\/2 
\z\/2 

V2 



dz. 



+ 



[£>,(r) + 2)9]V: 



■ exp 



2(£>,(r) + 2/3) 



[/3 + J D 2 (oo)- J D,(r)/2]V2 
Particular case of eq. B3 is the second moment 



exp 



(v+ - v ga i,+f 



f3 + D z (oo)-D z (r)/2 



(B3) 



{p s (X 1 ,v 1 ] 



dz 1 / dz 2 (p(X 1 ,z 1 )p(X 1 ,z 2 )) {</>vipti, ZiWviQZ-u Z2)) 



2vr /_ 



= — dz 



X 



S-\z\/2 

z\/2 

V2 



dz 4 



*(M) 



[D,(\z\) + 20\V* 



exp 



J gal 



2{D z {\z\) + 2(3) 



exp 



(vi - v gal ,+f 



f3 + D z (oo)-D z (\z\)/2 



(B4) 



[f3 + D z (oo)-D z (\z\)/2]^ 

Here we have used our standard notation v = V\ — v 2 , v + = (vi + v 2 )/2, z — z\ — z 2 , 
z + = (z! + z 2 )/2, r = xi - x 2 , Vg a i = u fla j(xi) - v gal (x 2 ) and ^ aij+ = [*v( x i) + v gai{^2)}/2. 

We observe that the density in the velocity space is statistically inhomogeneous which is 
reflected first of all in a residual dependence of the quantities in eqs. B2-B4 on the absolute 
velocity V\ or v + and not only on velocity difference v. In particular, eq B2 describes the 
mean velocity profile of the of the density, which is not, in general, uniform. 

There can also be dependence of the statistics on sky coordinates X, if the regular flow 
pattern v ga i(xL) is complex. We shall not consider this possibility, restricting our attention 
either to v ga i described by linear shearing pattern or the case where regular flow can be 
neglected altogether as is the case for isolated clouds. 

Expressions eqs. B2-B4 are quite complex. They are significantly simplified in several 
astrophysically important cases. First of all, when one considers small enough scales so that 
turbulent velocities are smaller than Z^(oo) 1 / 2 , all exponential terms which contain D z {po) 
can be taken equal to unity and we have 



(psiX^vi}) = const 
(psiXi^i) 2 ) = const 

(p a (Xi,Ui)p 8 (X 2 ,U2)) oc 



exp 



Vnalj 



2(D z (r) + 2p) 



(B5) 
(B6) 

B7) 



SJ [D,(t) + 20\V* 

This equations are applicable when the velocity line is extended compared with the turbulent 
scales under study which is the case, for example, for the 21cm emission of interstellar HI. 
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Henthforth we shall not follow precise values of the proportionality factors, which depend 
on the extend of the line-of-sight integration S and the details of the integral cut-off. For 
the analysis of the scaling laws at small separations R the boundary effects can be neglected 
and the range of integration extended to infinity. 

In another limit, when dealing with lines from individual clouds, one can consider regular 
flow to be absent. In this case 



(p s (Xi,wi)) = 
(p a (X 1 ,v 1 )p 8 (X 2 ,v 2 )) oc 



pS 



^[D z (oo) + 2f3} 1/2 
P 2 S 



exp 



D z (oo) + 2f3 



x 



(/9 a (Xi,Ui) 2 ) OC 



[D g (oo)+0\V* 
s 

dz 

-s 

p 2 S 2 

[D z (oo) + 0\V* 



exp 



e(r)/p 2 



(B8) 
(B9) 



[D,(r) + 20\V* 



exp 



2(D z (r) + 2(3) 



exp 



£> 2 (oo)+/3 



(BIO) 



CO. 



Again, for small separations R, the integration can be extended to an infinite range, S 

An image of an individual cloud in PPV space has an extension ~ D 1 ^ 2 (S). This 
stems from the fact that for an individual cloud the motions on the large scale determine 
the velocity dispersion. The assumed distribution, e.g. the Gaussian distribution, of large 
scale cloud velocities is present for an ensemble of clouds, while for an individual cloud 
under study its extension in the velocity space is a fixed realization from this distribution. 
Typically, for an individual cloud, the mean density is fixed. As the result the terms like 
exp[— v 2 /2{(5 + D(S))] in Eqs. (B8-B10) represent numerical factors and will be omitted in 
our further discussion. 



C. Correlation and structure functions in PPV 

The study of the 3D correlation function in PPV space has been performed in LP00, 
with an emphasis on the case of galactic turbulence characterized by the presence of the 
coherent flow due to galactic rotation. These results follow from eq. (B7) which can be 
written in the form equivalent to eq. (CI) of LP00 if temperature broadening is taken to be 
small (or considered later as a part of the velocity slice width): 



&(R,u) = £ v R,v) + Z P (R,v) 

f°° P 2 
£ v (R,v) oc J dz - exp 



(CI) 



Vgol) 



2D z (r) 
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i P (R,v) oc 



dz 



f(r) 



£>,(r)V2 



cxp 



2D z (r) 



In this Appendix we present complimentary set of results for the PPV structure func- 
tions, suitable for description of the gas confined to isolated clouds. Our starting point is 
the eq. (B9). We shall be interested only in the small scale r/5 1 < 1 behaviour, for which 
we use simplified expression 



£„(R,u) oc 



i p {R,v) oc 



P 2 S 



p 2 S 



L 



" i 

dz n 7 u/2 ex P 
s D z (r)y 2 

dz n 7 m/2 exp 
-s D z (v) 1 ' 2 



2D z (r) 

v 2 
~2Djv) 



(C2) 
(C3) 



D z (sy/ 2 

where boundary inhomogeneous effects are neglected. The quantity p s ps pS/D z (S) 1 ^ 2 is the 
approximate mean density in PPV space. At small amplitudes of turbulence velocity one 
should restore the thermal term in this expression, p s pS/(D z (S) + 2/3) 1 ^ 2 

Observations of interstellar turbulence testify that the emissivity obeys power-laws (see 
Armstrong, Rickett & Spangler 1995, Lazarian 1999, CLV02). This makes it natural to 
assume that the underlying velocity and density statistics is also power-law. Power laws for 
the statistics of turbulence are also expected from theoretical considerations (see CLV02). In 
what follows we shall use normalized variables, which for a cloud of extend S would amount 
to using spatial coordinates normalized by S and the velocities normalized by the velocity 
dispersion at the cloud scale CS m . For simplicity, we shall also write D z (r) = Cr m , ignoring 
non-critical additional angular dependence exhibited by pure solenoidal or potential flows. 
Asymptotics for d v (R,0), d p (R,0), R/S -> 0, as well as d v (0,v), d p (0,v), v 2 /CS m -> 0, 
follow from analysis of the integral (C3). Namely 

d v (R,0) = 2[&(0,0)-£ t ,( J R,0)]cx 



p 2 S 2 
D Z (S) 



oc 



oc 



dz 



z m ' 2 [(R/S) 2 + z 2 ] m / A 



oc 



d p (R,0) = 



oc 



p 2 S 2 y/HXm/4-1/2) 
D Z (S) r(m/4) 

2[O(0,0)-O( J R,0)]oc 

p 2 S 2 (r /Sy 



(R/S) 



l-m/2 



(C4) 



oc 



D Z (S) 
p 2 S 2 (r /Sy 



dz 



1 



1 



z 7+m/2 

r( 7 /2 + m/ 



[(R/S) 2 + 2 2]7/2+ m /4 

1 - 1/2) fR> 1 " W2 
S 



oc 



D Z (S) r( 7 /2 + m/4) \SJ (C5) 

The latter estimation is valid for 7 + m/2 < 1, i.e for sufficiently steep density spectra in real 
space. Kolmogorov density spectrum 7 = —2/3 belongs to this class. For shallow spectra 



-42- 



which have positive 7, one should use the correlation function £ p in place of d pi but one finds 
(see LPOO) that £ p (R, 0) still obeys the scaling (C5). 

Somewhat more complex is the velocity dependence of the structure functions. For 
simplicity, we shall describe all the steps in detail for the velocity term. 



d v (0,v) = 2[&(0,0)-£»]<x 
p 2 S 2 i- 1 



oc - ; — / dz 



oc 



D Z (S) J_ x W 2 
p 2 S 2 1 



1 — exp 



v 2 



2z r ' 

[l- P 2~ p v 2p T(-p,v 2 /2)} (C6) 



D Z (S) mp 

where p — 1/m — 1/2 > and V is incomplete gamma-function. For small argument, the 
series expansion of the incomplete gamma-function is 

(op op-1 \ 
- + v 2 + 0(v A ) (C7) 
p 1 -p J 

Substituting (C7)into (C6) we see that the most divergent term in the T[—p,v 2 /2] cancels 
out, while all the residual terms are convergent. The leading behaviour for the velocity 
structure function is d v (0,v) ~ r(— p)v 2p for p < 1 which changes to d v (0,v) ~ 4(^-1) ^ 2 f° r 
p > 1. Thus the structure function scaling depends on m for m > 2/3, d v (0,v) ~ v i 2 - m )/ m 
and is fixed to d v (0,v) ~ v 2 at m < 2/3. 

4(0,.) oc ^r(---)y , m >2/3 (C8) 

*m k S^ii • m<2/3 (C9) 

Notably, Kolmogorov value m = 2/3 is the boundary case here. One should note that 
representation (C7) is incorrect for integer p, i.e for m = 2/3,2/5,2/7. . .. Series expansion 
of the incomplete gamma function is irregular in these cases, which manifests itself in the 
appearance of logarithmic terms. For example 

T[-l,v 2 /2} ~ 2v~ 2 + logv 2 + O(v ), m = 2/3 (CIO) 
T[-2,v 2 /2] ~ 2v~ A -2v~ 2 -0(}og{v)), m = 2/5 (Cll) 

Thus, quadratic scaling of the structure function with velocity is preserved even in these 
cases, except that for the Kolmogorov turbulence, logarithmic correction enters the leading 
behaviour. 

As in the case of the radial dependence (C4), the velocity term can be considered 
a special case 7 = of the density term. Indeed d p (0, v) is obtained analogous to just 
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completed analysis by substituting p — 1/m — 7/m — 1/2 and restoring (ro/S*) 7 amplitude 
factor. 

1-7 1 

* |J(F r G- i ^ 2 )(Di)"'"- ™>2/3d-7)(C12) 

K H) 2d--, 1 ) -3m ^ • m < ^ " 1) < C13 » 

Notably, the range of m where velocity scaling is sensitive to m decreases for the long-range 
dominated density spectra, 7 < 0. For example, for Kolmogorov density spectra 7 = —2/3, 
we should use the second scaling solution for all m < 10/9. Thus, the quadratic scaling 
solution (C13) is the one of predominant interest in the study of such density inhomogeneities. 
On the other hand, for the short-wave dominated density spectra 7 > the m sensitivity is 
retained over much wider range and one should use expression (C12). 17 Exact transitional 
value ofm = 2/3(l — 7) points again to the presence of the slow varying logarithmic term 
in the asymptotic scaling. 

One may think that there is a contradiction of our present results with those in LP00. 
Indeed, spectra along the v-direction obtained in LP00 were found to reflect the spectral 
index over the whole range of spectral indexes studied. At the same time, structure functions 
of fluctuations measured along the velocity coordinate do not have this property. This 
contradiction resolves in a very simple way. The structure functions cannot grow faster than 
v 2 . As the result they do saturate at this value. At the same time spectra do not have such a 
limitations and properly reflect very steep spectra that arise. Therefore to identify the power 
index of fluctuations along the velocity coordinate it is necessary to measure spectra along 
the velocity axis. This provides a new way 18 to study turbulence that is the subject of a 
current research in collaboration with Alexey Chepurnov that we shall elaborate elsewhere. 
Here we shall only note that the steepness of the spectrum prevents a useful generalization 
of the Spectral Correlation Function approach (see Rosolowsky et al. 1999) to the velocity 
direction (compare to a discussion in Lazarian 1999). 

17 Howcvcr, for 7 + m/2 > 1 our analysis should be replaced by the PPV correlation function formalism. 

18 Note, that it is not necessary to have good spatial resolution to get turbulence if spectra along the velocity 
coordinate are taken. Therefore velocity turbulence studies may become an important part of extragalactic 
research. The corresponding asymptotics can be found in LP00, Table 3. 



-44- 



D. Optically thin case: structure function approach 

In this Appendix we revisit the problem of the intergalactic turbulence in the case of 
optically thin lines, which was the focus of LPOO. With the correlation function in velocity 
space given by (18) we are now able to directly derive the asymptotic expressions for the cor- 
relation (or structure) functions of the intensity, bypassing rather cumbersome calculations 
of the 3D power spectrum in the velocity space. 

The intensity of the emission in an optically thin line is proportional to the density of 
the atoms in the velocity space 

I(X u v 1 ) = ep a (X 1 ,v 1 ) (Dl) 



The mean line profile follows from eq. (C3) 



(/(Xi,t;i)> = 



ep 



(2tt[/3 + D 2 (oo)]) 



1/2 



dz exp 



2[P + D z (oo)} 



(D2) 



Total intensity in a velocity channel of width SV centered at velocity V and described 
by the shape function Wc(v — V), such that / dvWc(v — V) = 5V is 



I c (Xi) = e J dv, W c ( Vl - V) p s (X 1 ,v 1 ) 



The mean of the channel-integrated intensity is 



ac(xo) 



ep 



(2vr [P + D z (oo)]) 1/2 
For the correlation function, 



l dz I 



dvWc(v — V) exp 



2[/? + J D 2 (oo)] 



(D3) 



(D4) 



6(R) = (/c(Xi)/ c (X 2 )> = e 2 J dv 1 W c (v 1 ) J dv 2 W c (v 2 ){p s (X 1 ,v 1 )p s (X 2 ,v 2 )) , (D5) 

we shall use expression (B7) applicable to the case of intergalactic turbulence with extended 
gas distribution and coherent galactic flow with a linearized shear: 



/c 



dz 



l + |(r) [D z (r) + 2(3\ 



-1/2 



dv W{v ) exp 



(v - r i zf 



2(D z (r) + 2(3) 



(D6) 



Here two velocity channel windows are combined into W{v) = f dv + Wc(v + /2—v)Wc(v + /2+ 
v). 
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Equation (D6) describes the effect of the line-of-sight projection of 3D turbulent distri- 
butions of emitters. Velocity kernel both modifies the projection of the density correlation 
function (£(r) term in the first bracket), which now exhibits different scaling than just a 
column density correlation function would, and also introduces a new phenomena - even 
uniform spatial distribution of emitters following incompressible turbulent flow will be cause 
with fluctuating intensity if observations are carried out in sufficiently narrow velocity chan- 
nel. The pure velocity induced fluctuations are defined by the first, namely, unity term in 
the first bracket. Further one, we shall discuss two effects separately, using obvious split of 
the full correlation function into two terms £/(R) = £y(R) + £p(R)- 

As expected, the intensity integrated over the whole line contains no information about 
the velocity field. Indeed, integration with a constant window W{y) over all velocity range 
in eq. (D6) gives 

^ (R) ~2^ J J*^^ Rl ~" (D7) 

while is reduced to a constant which describes uniform mean intensity distribution and 
can be dropped. 

In LPOO we presented a detailed discussion of the velocity influence on the projected 
intensity statistics in optically thin lines. Here we shall reproduce the scaling solutions at 
small sky separations R which was shown in LPOO can be found for observations in thin 
velocity channels, using correlation function formalism of eq. (D6) as the starting point. 



For the very narrow channel W(v) = 5(v) we have 

1 f°° r i 
f (R) ~ _ J dz [l + |(r)J [D z (v) + 2/T 1/2 exp 



2(L> 2 (r)+2/3) 



(D8) 



The asymptotic behaviour of the density term at the small scales R <C (C 2 /) 2 - m is obtained 
st r aight forwar dly : 

-i poo -i poo 

UK)~-J azl(r)[D z (T) + 2P]-^~-J dz t{v)/D z {r)-^~ R^~^ . (D9) 

The latter expansion, valid when R is still large enough for the turbulent dispersion to exceed 
the thermal one D(R) > (3, is the one of the thin slice regime. Importantly, the slope is 
modified by the term m, reflecting the statistical properties of the turbulent velocity. 

Asymptotic analysis of the pure velocity contribution 



£,(R) ~ ^ / <l: [ D -( r )' 2, l" 1/2 ^ 



2(D 2 (r) + 2 / 3) 



(DIO) 
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case f(R) P(K) 

totcifc slice R 1 ^ K~ 3+ i 

thin slice R l-l-m/2 x - 3+7 +m/2 

tain s/ice, ve/oc% R}~™l 2 K~ 3+m / 2 



is more delicate, since the integral dz [D z (r) + 2f3]~ 1 / 2 does not converge and the expo- 
nential factor cannot be just set to unity in the leading order. The following manipulations 
lead to the correct result: 

Using D(r) = C(R 2 + z 2 ) m ^ 2 , omitting thermal term for brevity and introducing new 
variable y = \ m l 2 ~ l zj{R 2 + z 2 ) m ^, A = (C 1 / 2 /)2/2-m we obtain 



&(R) 



2tt 



i? 2 + ^ 2 



-y 2 /2 



R 2 + z 2 (l -m/2) 

where we have left z as the implicit function of y. Now we integration by parts 



(Dll) 



&(R) ~ £ ^ erf 



y 

V2. 



R 2 + z 2 



R 2 + z 2 (l - m/2) 



dy erf 



1/ 

v/2. 



R 2 + z 2 



Ay R 2 + z 2 (l - m/2) 

(D12) 

The first term is finite, while in the second one we replace the integration variable back to z — 
z/R and expand the error function in power series in y = (R/ \) l ~ m ^ 2 z/{l + z 2 ) m ^ keeping the 
dominant first linear term. The remaining integral is convergent and the scale dependent 
term in the correlation function can be expressed through Gamma and Hypergeometric 
functions 

t (j3 \ 2 _ f ff /ui-m/2 T[m/4-3/2] /3 5 m m\ 

UR) WX) 2 3 / 2 T[m/4] 2Fl U' ' 2 " T' 1 " T ) ■ (D13) 

Actually, only m/(2 — m) part of the second moment £ v (0) = 2/(2 — m) = 1 + m/(2 — 
m) describes the variance of the fluctuations of intensity, while the extra unity (in our 
normalization) comes from a nonzero mean value of /(X). The structure function, insensitive 
to the mean value of the field, is 



DJR) ~ 



(R/X) 



l-m/2 



r[m/4-3/2] 
2V2r[m/4] 



5 



m 



m 



F 1 - 2,---,!-- 



(D14) 



The Table D summarizes these small i? asymptotics of correlation functions, and the 
corresponding behaviour of the power spectra 
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The correlation function formalism, developed in this Appendix, also makes self-evident 
the distinction between thin and thick channels, introduced in LPOO. Direct inspection of 
eq. (D6) shows that one can recover thin slice behaviour if the width of the channel window 
5V combined with thermal broadening does not exceed the turbulent velocity amplitude 
at the scale of consideration 5V 2 + 2/3 <C D(R) while if opposite is true the slice is thick. 
Naturally, this criterion coincides with that derived in LPOO. 



E. Gaussian fluctuations: General case 

While in section 3 we attempted to be as general as possible while dealing with corre- 
lations at small angular separations, here we deal with arbitrary angular separations. The 
price that we pay for that is that we have to make an assumption about the statistics from 
the very beginning. Namely, we assume that fluctuations are Gaussian with the probability 
function: 

1 ( 1 

C = A 2 -e s , A={(6p 3 y), (El) 
A is the dispersion and the information about the field correlation is contained in £ s . 

The exponent in (26) can be rewritten in terms of fluctuations exp(ap s ) = exp(ap s ) exp(aSp s 
Then eq. (32) can be rewritten: 



P(R) = ^ J dvxWivt) J dv 2 W(v 2 ) 



x 



x e 



-a(pn+pi 2 ) |^ e -ct(5pii+5pi2) _|_ g -a(<5p22+5p2i) _ g -a(5p2i+5pi2) _ g-a(<5pn+<5p22)^ J(E2) 



To evaluate the averages in (29) with probability distribution (El) we observe 19 that 



(exp [-a(5pu + Spjj)]) = exp 



— ((Spii + Spjj) 2 ) 



(E3) 



where quantities 5pu and 5pjj are measured at points Xi,fj and Xj,t>j of the data cube. 
Assuming homogeneity in the velocity distribution we get correlations that depend only 



19 The procedure of averaging with distribution given by eq.(El) is straightforward (see LPOO) and amounts 
to taking Gaussian integrals over yi and 7/2- 
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on the velocity difference v 20 between the points and the separation between the points: 

P(R) oc ^ J dvW 2 (v) (exp [a%(0,v)] - exp [a 2 £ s (R, *;)]) , (E4) 

The expansion of the exponent for small separations R results in eq. (37) taking into account 
that 2[&(0, v) - &(R, «)] = 4(R, u) - 4(0, v). 
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